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Abstract 


The central object of study of this thesis is inverse mean curvature vector flow of 
two-dimensional surfaces in four-dimensional spacetimes. Being a system of forward- 
backward parabolic PDEs, inverse mean curvature vector flow equation lacks a gen¬ 
eral existence theory. Our main contribution is proving that there exist inhnitely 
many spacetimes, not necessarily spherically symmetric or static, that admit smooth 
global solutions to inverse mean curvature vector flow. Prior to our work, such solu¬ 
tions were only known in spherically symmetric and static spacetimes. The technique 
used in this thesis might be important to prove the Spacetime Penrose Conjecture, 
which remains open today. 

Given a spacetime (A^, g) and a spacelike hypersurface M. For any closed surface 
S embedded in M satisfying some natural conditions, one can “steer” the spacetime 
metric g such that the mean curvature vector held of S becomes tangential to M 
while keeping the induced metric on M. This can be used to construct more examples 
of smooth solutions to inverse mean curvature vector how from smooth solutions to 
inverse mean curvature how in a spacelike hypersurface. 
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1 


Introduction 


1.1 Motivation: Mass in General Relativity 

General relativity is the study of large scale structures of the universe. One funda¬ 
mental object in general relativity is the notion of mass. Pointwise energy density 
and total mass of a spacetime are both well-dehned in general relativity. However, 
the local mass of a given region in a spacetime (called quasi-local mass), as well as 
the relationship between local mass and pointwise energy density and total mass of 
the spacetime are still not very well understood. 

Despite of many attempts in dehning the quasi-local mass (e.g. [2, 3, 5, 15, 16, 
24, 47]), none of the proposed functionals satisfy all the desired properties. One 
such natural property is that the total mass of the spacetime should be bounded 
from below by the mass of a region in it, assuming some positivity condition on the 
pointwise energy density (e.g. dominant energy condition). 

Given a spacetime {N‘^,g) and a complete asymptotically flat spacelike hypersur¬ 
face (also called a slice) with the induced Riemannian metric g. Let k be the 
second fundamental form of M. The triple {M^,g, k) is called a Cauchy data of this 
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hypersurface (see Figure 1.1 below). 


■A 

p 

V 

k) 

Figure 1.1; Slice in a spacetime with Cauchy data. 

There is a well-dehned quantity called the ADM mass (dehned by R. Arnowitt, 
S. Deser and C. Misner in [1]) that measures the total mass of this hypersurface. 
Suppose M has a compact outermost minimizing surface S. Physically, S can be 
viewed as the apparent horizon of blackholes. 

In the case that M is totally geodesic, i.e. fc = 0, then the pointwise energy density 
equals the scalar curvature of M. In this case, the Riemannian Penrose Inequality 
states that: 

Theorem 1.1 (Riemannian Penrose Inequality). Let M and S be given as above. If 
the scalar curvature of (M, g) is non-negative, then its ADM mass is greater than or 
egual to , where |S| is the total area ofH (see Figure 1.2). 

Penrose [37] hrst conjectured this inequality in 1973, and he gave a heuristic proof 
based on physical considerations, explained as follows. It turns out that the lower 
bound -y/lSl/lbTT in the Riemannian Penrose Inequality equals the Hawking mass, 
which is a quasi-local mass functional proposed by Hawking [24], of the minimal 
surface S. This can be viewed as the mass of the blackholes inside S. Thus, the 
Riemannian Penrose Inequality states that the total mass of M should be at least the 
mass contributed by the blackholes, assuming that the energy density (which equals 
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Figure 1.2: Totally geodesic spacelike hypersurface with compact outermost mini¬ 
mal surface S. 

the scalar curvature in the case of a totally geodesic hypersurface) is non-negative 
everywhere. 

Geroch [22], Jang and Wald [26] hrst discovered a monotone property of Hawking 
mass of surfaces under smooth inverse mean curvature flow. Based on this, Huisken 
and Ilmanen [25] gave a proof of this inequality in the case of a single blackhole (i.e. 
S is connected). In the same year, Bray [5] proved the full Riemannian Penrose 
Inequality using a different technique. 

In the case of no blackholes, the Riemannian Penrose Inequality is also known as 
the Riemannian Positive Mass Theorem: 

Theorem 1.2 (Riemannian Positive Mass Theorem). Given a complete asymptot¬ 
ically flat Riemannian manifold with non-negative scalar curvature. The 

ADM mass of M is non-negative. 

In 1979, Schoen and Yau [42] proved this result using a variational method. In 
the same year, they [41] generalized this result to Riemannian manifolds of dimen¬ 
sion less than eight. In 1981, they [43] removed the assumption that M is totally 
geodesic and proved the Riemannian Positive Mass Theorem for an arbitrary space- 
like hypersurface in a spacetime that satishes the dominant energy condition. 
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1.2 Inverse Mean Curvature Vector Flow 


So far all the discussions assume that our spacetime has a hypersurface with zero 
second fundamental form, in which case inverse mean curvature flow naturally bridges 
the Hawking mass of an apparent horizon of blackholes and the ADM mass of the 
hypersurface. However, a spacetime in general does not admit such a totally geodesic 
hypersurface. This is because that the second fundamental form of a hypersurface 
in a spacetime has six components, but the hypersurface only has one degree of 
freedom. Thus it is not generic to have all six components vanish. Therefore it is 
desirable to obtain a similar bound on the total mass of the spacetime by the mass of 
the blackholes without this assumption. This leads to the general Spacetime Penrose 
Conjecture, which is still open today. 

A viable candidate for proving this conjecture is the codimension-two analogue 
of inverse mean curvature flow, called the inverse mean curvature vector flow. How¬ 
ever, there are two major problems with this flow. First, unlike the inverse mean 
curvature flow, which is a forward parabolic PDF, inverse mean curvature vector 
flow is a system of forward-backward parabolic PDFs: forward parabolic in space- 
like directions, and backward parabolic in timelike directions (see [25]). Backward 
parabolic equations lack a general existence theory. For instance, the reverse heat 
flow is backward-parabolic. Given many initial conditions, the reverse heat flow 
would develop singularities instantaneously. However, the reverse heat flow would 
exist for time f > 0 if we hrst perform the heat flow for time t and then start flowing 
backwards. 

Second, for some initial surfaces, even the inverse mean curvature vector flow 
exist, their Hawking mass still won’t give us a lower bound on the total mass of the 
spacetime as in the inverse mean curvature flow case simply because the former is 
too large. To illustrate this, take a t = constant slice in the Minkowski spacetime. 
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The round sphere in that slice has zero Hawking mass. Spacial perturbations will 
decrease the Hawking mass making it negative, whereas timelike perturbations will 
increase the Hawking mass making it positive. During inverse mean curvature vector 
flow, the spacial “wiggles” will smooth out due to the parabolic nature of the flow. 
However, timelike “wiggles” will get amplified since the flow is reverse parabolic in 
the timelike directions. With these surfaces with positive Hawking mass, inverse 
mean curvature vector flow will not provide a lower bound for the ADM mass of 
Minkowski space, which is zero. 

However, these two problems seem to solve each other because they are both 
suggesting that solutions to inverse mean curvature vector flow exist only when the 
“right” initial surface is given. The important question is then: Given a space- 
time. Do such “right” initial surfaces always exist? The answer is affirmative if the 
spacetime is spherically symmetric or static. 

Inverse mean curvatnre vector flow of snrfaces in spherically symmetric space- 
times was first studied by E. Malec, and N. OMurchadha [31]. They showed that 
inverse mean curvature vector flow of spherically symmetric spheres exist for all time. 
Intuitively, the spherical symmetries prevent timelike “wiggles” to occur. Moreover, 
the Hawking mass is monotonically non-decreasing under inverse mean curvature 
vector flow of spacelike surfaces with spacelike mean curvature vectors, assuming the 
spacetime satisfies the dominant energy condition. 

Later Frauendiener [21] showed that, in an arbitrary spacetime that satisfies the 
dominant energy condition, if smooth inverse mean curvature vector flow exists, then 
the Hawking mass is monotonically non-decreasing. In 2004, H. Bray, S. Hayward, 
M. Mars and W. Simom [8] showed that we can in fact flow along a one-parameter 
family of directions and the Hawking mass is still monotone. 

In spherically symmetric spacetimes, the “right” initial snrfaces for inverse mean 
curvature vector flow are spherically symmetric spheres. What about spacetimes 
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that are not necessarily spherically symmetric? 

Bray and Ye Li were trying to develop a general existence theory for inverse mean 
curvature vector flow back in 2009, and one of their intuitions was that if one can 
somehow control the flow of the surfaces so that they stays tangential to a spacelike 
slice, then the flow might not develop singularities. In fact it has been shown that: 

Proposition 1.3 ([8]). The family of closed embedded spacelike surfaces is a 
solution to the smooth inverse mean curvature vector flow with spacelike inverse mean 
curvature vector everywhere on the surfaces if and only if there exists a spacelike 
hypersurface c N, such that the mean curvature vector is tangential to M 
at all {x, s), and {Ss} is a solution to the smooth inverse mean curvature flow in M. 

Following the intuition, we prove the main theorem in Chapter 4: 

Theorem 1.4 (Main Theorem). There exist infinitely many non-spherically sym¬ 
metric, non-static spacetimes that admit inverse mean curvature vector flow coordi¬ 
nate charts. Given such a spacetime U with an inverse mean curvature vector flow 
coordinate chart (t, r, 9, f) and the constructed spacetime metric g. The coordinate 
spheres St^r contained in each t = constant slice, when reparameterized by r^ = e^, 
are smooth global solutions to the inverse mean curvature vector flow equation. 

This theorem is restated and proved in Theorem 4.4. The proof is based on 
explicit constructions of spacetime metrics that admit inverse mean curvature vector 
flow coordinate charts, dehned in Chapter 4. Theorem 1.4 seems to suggest that 
spacetimes that admit smooth solutions to inverse mean curvature vector flow exist 
generically. However, this general problem of hnd solutions to inverse mean curvature 
vector flow (i.e. the “right” initial surface) in arbitrary spacetimes is still open. 

There also exists a coordinate-free analogue of Theorem 1.4: 
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Theorem 1.5. Given a spacetime a spacelike hypersurface and a closed 

embedded surface c M. Suppose S is area expanding (defined in (4-3.5)), 

then there exists a unique smooth steering parameter Q = Qs ^ C°°(iV), such that 
in the steered spacetime metric gq (defined in (4-3.1)), H^, is tangential to M every¬ 
where on S. 

This can be used to generate more examples of solutions to inverse mean curvature 
vector flow. Consider a smooth solution to inverse mean curvature flow in a spacelike 
hypersurface M. One can then smoothly adjust the spacetime metric along the flow 
such that the mean curvature vector of each flow surface becomes tangential to M. 
These steered surfaces are then solutions to inverse mean curvature vector flow, since 
the area expanding condition is already satisfied, and now the mean curvature vectors 
are tangential to a spacelike hypersurface (see Proposition 1.3). 

1.3 Uniformly Area Expanding Straight Out Flows and Time Flat 
Surfaces 

Inverse mean curvature vector flow is a type of flow that has bad existence theory, 
but very good properties: the Hawking mass is monotone under smooth inverse mean 
curvature vector flow. There is another flow studied in Chapter 5, called uniformly 
area expanding straight out flow (or simply straight out flow), that has solutions with 
a wide class of initial surfaces. In that chapter, we try to construct spacetimes that 
admit a coordinate chart in which straight out flow of coordinate spheres exists for 
all time. Partial results have been obtained while complete understanding of this 
problem is still work in progress: 

Proposition 1.6. Suppose a spacetime {N^,g) admits a coordinate chart {f,r, 0,0} 
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such that the coordinate representation of g is 


t 



r 

d 


9 0 



9 = 

6 


d 0 0 

e 0 a c 


(j)\ f 0 c b J 


(1.3.1) 


with ab — = r^siv?9. Then := is straight out if and only if d satisfies a 

second order elliptic PDE in d: Ag^d + G{d,d') = 0, where G is given by (5.2.30) 
and (5.2.31). 


We conjecture that solutions to the above elliptic PDE always exist. 

In addition to general existence, another reason for studying the straight out flow 
is that the Hawking mass is also monotonically non-decreasing under such flow (e.g. 
see [9], [10]) if the spacetime also satishes the dominant energy condition. 

The general existence of straight out flows can serve as a disadvantage since 
we can even flow surfaces with positive Hawking mass, too large to be used as a 
lower bound of the total mass of some spacetimes, in straight out directions. To 
see this, again consider the Minkowski spacetime. All surfaces that are contained in 
a spacelike plane have non-positive Hawking mass. Thus, non-planer surfaces have 
positive Hawking mass. For such surfaces, inverse mean curvature vector flow would 
not work since there are time “wiggles”. However, those surfaces can still flow in 
straight out directions. Since the total mass of Minkowski space is zero, having a 
surface with positive Hawking mass is not going to give us a lower bound for the 
total mass since the Hawking mass is monotone. 

Inverse mean curvature vector flow and uniformly area expanding straight out 
flow are two special cases of uniformly area expanding flows: orthogonal flows such 
that the rate of change of the area form of each flow surface equals the area form 



itself. The Hawking mass is not necessarily monotone under general uniformly area 
expanding flows. H. Bray, J. Jauregui and M. Mars very recently ([9], [10]) obtained 
a variational formula of the Hawking mass under general uniformly area expanding 
flows, which consists of four major terms (see [10]). The hrst three terms are non¬ 
negative if the spacetime satishes the dominant energy condition. The fourth term 
is an integral term with integrant a function of the spacetime multiplied by the 
divergence of the connection one-form associated with their mean curvature vector 
of the flow surfaces. Thus, if the connection one-form is divergence free, then the 
fourth term vanishes and the Hawking mass is monotone. Surfaces with divergence 
free connection one-form associated with the mean curvature vectors are called time- 
flat (dehned in [9], [10]). While the conditions on inverse mean curvature vector flow 
coordinate chart can be viewed as a global “flatness” condition on the surface, the 
time-flat condition is a local “flatness” condition. 

The organization of this thesis is given as follows. In Chapter 2, we study the 
monotonicity of Hawking mass under smooth inverse mean curvature flow. In Chap¬ 
ter 3, we study inverse mean curvature vector flow in spherically symmetric space- 
times. The notations used in this thesis are also introduced in that chapter. In 
Chapter 4, we prove the main theorems 1.4 and 1.5. In Chapter 5, we study uni¬ 
formly area expanding straight out flows, and prove Proposition 1.6. Finally in 
Chapter 6, some open problems and future works are discussed. 
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2 


Huisken-Ilmanen Inverse Mean Curvature Flow and 

Monotonicity of Hawking Mass 


In this chapter we study inverse mean eurvature flow of a closed embedded surface 
in an asymptotically flat Riemannian manifold {M^,g). The motivation is the 
Riemannian Penrose Inequality: 

Theorem 2.1 (Riemannian Penrose Inequality). Let {M^,g) be a complete, asymp¬ 
totically flat Riemannian manifold with non-negative scalar curvature and a compact 
outermost minimal surface S of total area |S|, then 


rOADM ^ 


IGtt 


( 2 . 0 . 1 ) 


with equality if and only if {M^,g) is isometric to the Schwarzschild metric with mass 
m > 0.' 


m 


KAio}. + 




( 2 . 0 . 2 ) 


outside their respective outermost minimal surfaces. 
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S here can be viewed at the apparent horizon of blackholes. The lower bound for 
the ADM mass, has the physical interpretation as the mass of the blackholes. 

Penrose [37] first conjectured the Riemannian Penrose Inequality in 1973, and he 
gave a heuristic proof based on the physical considerations. In 2001, Huisken and 
Ilmanen [25] proved this inequality using inverse mean curvature flow in the case of a 
single blackhole. In the same year, H. Bray [5] proved this inequality using conformal 
flow of metrics that works for any number of blackholes. In 2009, H. Bray and D. 
Lee [13] generalized the inequality to all dimensions less than eight. In 2010, Lam 
[28] proved the Riemannian Penrose Inequality for graphs in all dimensions. In 2011, 
Schwartz [45] proved a volumetric version of the Penrose inequality for conformally 
flat manifolds. The general Spacetime Penrose Conjecture is still open today (see 
[33, 11, 12, 33] for more discussions of this conjecture). 

In Section 2.1, we dehne asymptotically flat manifolds, ADM mass and Hawking 
mass of closed surfaces. In Section 2.2, we study the Geroch, Jang-Wald monotonicity 
formula of Hawking mass under smooth inverse mean curvature flow. In Section 2.3, 
we briefly discuss Huisken and Ilmanen’s proof of the Riemannian Penrose Inequality 
using such flows. 

2.1 Asymptotically Flat Manifolds, ADM Mass and Hawking Mass 

Definition 2.1.1. An n-dimensional Riemannian manifold is called asymp¬ 

totically flat if it satisfies the following two conditions: 

(1) There exists a compact set K ^ M and a diffeomorphism 

^:E:= M\K —^ 

where Bi is the unit open ball in M"; and 


11 



(2) In the coordinate chart ■ ■ ,x"') on E induced by the above diffeomor- 

phism called an asymptotically flat coordinate chart, the metric components 
Qij and the scalar curvature R satisfy the following decay conditions at any 
point X E E, i,j,k,I = 1,2, - ■ ■ , n: 

(1) piflx) = diflx) + 0{\x\-P); 

(2) \x\\gij,k{x)\ + \x\‘^\gipki{x)\ = 0{\x\-P); 

(3) \R{x)\=0{\x\-<^), 

with some constants p > and q > n. Here gij^k and gtj^ki are coordinate 
derivatives. 

E is called an asymptotically flat end of M. An asymptotically flat manifold can 
have multiple asymptotically flat ends. 

Definition 2.1.2. Given an asymptotically flat Riemannian manifold {M"-,g) and 
asymptotically flat coordinate chart. The ADM mass of of M is: 

If” 

mADM{M, g)-.= Wm. — --- 'V. ihijA - gii,j)nMSr (2.1.1) 

r^oo2(n- l)UJn-l 

where Un-i is the volume of the {n— 1) -dimensional round sphere; Sr is the coordinate 
sphere of radius r; v is the outward unit normal along Sr,' and dSr is the volume form 
of Sr. 

The ADM mass of an asymptotically flat manifold was defined by Richard Arnowitt, 
Stanley Deser and Charles W. Misner [1]. They proved that the above dehnition is 
independent of the choice of asymptotically flat coordinate charts. Thus, the notion 
of the ADM mass is well-dehned. We sometimes simply write mADM{g) instead of 
nrADM{M, g) if the underlying manifold is clear. 
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In dimension 3, we have 


mADM{M^,g) 


lim - 

r^co IGtT 


dS^. 

*0 = 1 


( 2 , 1 , 2 ) 


Definition 2.1.3. Given a Riemannian manifold {M^,g) and a closed embedded 
surface with the induced metric. The Hawking mass ofT is defined to be: 


: = 



(2.1.3) 


where H^. is the scalar mean curvature ofH in M. 


Example 2.1.1 (Euclidean Space). M” with the standard Euclidean metric is an 
asymptotically flat manifold with zero ADM mass. 


Example 2.1.2 (Conformal Transformation of Metric). Given an asymptotically flat 

4 

manifold {M^,g). Gonsider a conformal transformation g = u^-^g of the metric g, 
with u e u > 0. By Equation (A.5.14) the scalar curvatures R and R of g 

and g respectively, are related by: 

R = u-'^ [ru - ^9^ ■ (2-1-4) 


If u and its coordinate derivatives satisfy the following decay conditions, i,j,k = 
1,2, •• • ,n: 

(1) u tends to 1 at oo; 

(2) = 0{\x\-P-^); 

(3) ujk = 0(|a;|"P"2); 

(4) AgU = 0{\x\-^) 
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in an asymptotically flat coordinate chart of for some constants p > and 

q > n, then is also asymptotically flat in that coordinate chart. Moreover, 


mADAiig) = mADMig) - lim 7 -- -^dSr, (2.1.5) 

r-^oo (n - 2)Un-l Jsr 

where ^ is the outward normal derivative of u along Sr- 

Example 2.1.3 (Schwarzschild Manifold). Combining Example 2.1.1 and 2.1.2, con¬ 
sider the following one-parameter family of conformal transformations o/(M'^\{0}, Sij), 
parameterized by a constant m > 0: 

(l + ■5.,) . (2.1.6) 

This is called the Schwarzschild manifold of dimension n and mass m. It is easy to 
verify that m = 1 + 2\:^-2 satisfies the desired decay condition to make the resulting 
metric asymptotically flat, m is called the mass because the ADM mass of this metric 
is exactly m. This can be seen quite easily for the three-dimensional Schwarzschild 

manifold: . By Equation (2.1.5): 


rOADM 



mADM{d) - lim — 
r^oo ZTT 



dSr 


= — lim — 

r^oo 271 




dSr 


1- 1 -2i 2 

= lim-r 47 rr 

r^oo 27T 2 


= m 


(2.1.7) 


We now study further the geometry of the three-dimensional Schwarzschild man¬ 
ifold. Let r := |a;|. Eirst, note that the Schwarzschild metric is symmetric under the 
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2 

mapping r i—> ^. Thus the Schwarzschild manifold has two ends, with the center of 
symmetry being ^ that is r = ff, which is a two-sphere. Recall that if g = u'^S 

for some positive function u, then the mean curvature of a sphere of radius r with 
respect to g is given by: 


1 /2 Adu\ 
\r u dr J 


( 2 . 1 . 8 ) 


Therefore at r = 'j, the mean curvature is zero. Hence, the sphere r = ^ is a mini¬ 
mal surface, which can be viewed as the apparent horizon of a blackhole. The region 
outside of the blackhole is called the exterior region of the Schwarzschild manifold: 



The exterior region is an asymptotically fiat end (see Figure 2.1) 




Figure 2.1: Exterior region of three-dimension Schwarzschild manifold with bound¬ 
ary the minimal sphere Figure courtesy of Mau-Kwong G. Lam. 


There exists an isometric embedding of the three-dimensional Schwarzschild man¬ 
ifold into such that 

w‘^ 

r = -1- 2m. 

8m 

The image of this embedding is a parabola (see Figure 2.2), and the minimal sphere 
Sm. gets mapped to the sphere S 2 m ^ 
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Figure 2.2: Isometric embedding of three-dimensional Schwarzschild manifold into 
Mh 


The area of Sir is then given by the Euclidean area of S 2 m'- 1*5'^| = 47r(2m)^ = 
167rm^. Therefore 


mniSiR) 



Combing this with the ADM mass (2.1.7), we see that 


(2.1.9) 


Proposition 2.2. The ADM mass of the three-dimensional Schwarzschild manifold 
^]R^\{0}, -I- eguals the Hawking mass of the minimal sphere Sir, which 

is exactly m. 

More generally, Huisken and Ilmanen [25] proved (see also [40]): 

Theorem 2.3. Given an asymptotically flat Riemannian manifold and an 

asymptotically flat coordinate chart. Then 


lim mniSr) = mADM^M), (2.1.10) 

r^oo 

where Sr is the coordinate sphere of radius r. 
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2.2 Geroch, Jang-Wald’s Approach and Their Monotonicity Formula 


Huisken and Ilmanen’s proof of the Riemannian Penrose Ineqnality is based on the 
monotonicity property of the Hawking mass under smooth inverse mean curvature 
flow, first discovered by Geroch and Jang-Wald. 

Definition 2.2.1 (Inverse Mean Curvature Flow). Given a Riemannian manifold 
and a closed embedded surface in M. A smooth inverse mean curvature 
flow of Tj in M is a smooth family of surfaces F : S x [0, T] —> M of Yi such that 
the following parabolic evolution equation is satisfied: 

«no.n (2.2.1) 

where Vt and Ht is the unit outward normal vector field and scalar mean curvature 
of Yt := F{Y,t), respectively. 

A family of closed embedded surfaces {St} in M is called a smooth solution to 
inverse mean curvature flow if they satisfy (2.2.1). Given such a family of surfaces 
{St}. The hrst variation of area formula (A.3.7) implies that 

4|St| = f Ht^dAj: = |St|. (2.2.2) 

at Jzt ait 

Therefore the area of St grows exponentially under inverse mean curvature flow. 

Example 2.2.1 (Inverse Mean Curvature Flow of Spheres). Consider a round sphere 
SrQ in with radius ro > 0, and flow this sphere out by inverse mean curvature flow. 
By (2.2.2), the flow surfaces are still round spheres, and the area grows exponentially. 
Thus {Rgi/a^o} is a solution to this flow for all time. 

Geroch [22], Jang and Wald [26] discovered the following nice connection between 
solutions to inverse mean curvature flow and the Hawking mass: 
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Theorem 2.4 (Geroch, Jang-Wald). Given with non-negative sealar eurva- 

ture. If a family of elosed embedded surfaees {St} is a smooth solution to inverse 
mean eurvature flow in M, then for all t > 0, 

^m^t(St) ^ 0, (2.2.3) 

i.e. the Hawking mass is monotonically non-decreasing. 

Proof. Let Ht and dAt be the scalar mean curvature and the volume form of St in 
M, respectively. 


d d 


IGtt 


IhTT Js, ' * 


Hf dA 
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dt \ M IGtt 


IGtt 
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dAt I + 


1 d 


IGtt V IGtt dt 
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Ht dAt 


1 1 


2 Vievr 


,St 1- 

167r' ' V 167r 


Ht dAi 


St 


(By Equation (2.2.2)) 


+ 


I 


levr V 167r Js, 
1 . 1 


2Htj^{Ht)dAt + Htj^{dAt) 


2 V levr 




167r 2 


Ihvr js 
1 


Ht dAt ) + 


1 


IGtt V IGtt 


fL 

' dt ' 


2Ht^(Ht)dAt + HtdAt 


IGtt 


JEt 


Hf dAt 


1 

IhTT .;st 


d 




dAj 


(2.2.4) 


By the hrst variation of mean curvature formula (A.4.3): 




(2.2.5) 


where Ric^^ is the Ricci curvature of M. Plug (2.2.5) into (2.2.4) we get: 
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d 

dt 





1 

2 


1 1 
2 1^ 



H 


dAt 


IGtt 


2Ht 


St 



Ast -TT 


H, 


H, 


H, 


+ H^dAt 


1^1 

167rl2 ' IGttJs, 


1 1 
- + 


+ 2Ric"(!/,!/) + 2||II,|p - 


dAt 
( 2 ^ 2 . 6 ) 


We now compute the first three integral terms in the above. By integration by 
parts, we get: 

f 2H,M, (A) dA, - ( -2(y^,H,V^A)dA,= ( (2.2.7) 

Now by the Gauss equation (see e.g. [27]), we have 

Ric""(z/, z/) = - ^||Ih||2, (2.2.8) 

where is the Gauss curvature of E*. 

Next let \i(t) and A 2 (t) be the principal curvatures of E^, then 


Ht = \t{t) + \2{t), ||Ih||2 = Ai(t)VA2(t)'. 


Therefore 


- -Ht - \i{tY + \2{tY - 


[Ai(t) + A2(t)]2 [Ai(t)-A2(t)]' 


(2.2.9) 
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Plug them back into (2.2.6), we get: 




St L 


2i/,Aa I i) + 2mc’-<{y,y) + 2||II,||'" - 


E.lfl 1 


- + 


I 


16. I 2 ■ 16. I.. ^ ^ ^ 


> 


> 


Silfl 


167r\2 IGtt 

vit-it-".) 


r _ 2 A-& + 

Js, 2 


(i?" ^ 0) 


^ 0 


( 2 . 2 . 10 ) 


where the last inequality follows from the Gauss-Bonnet formula: 



K^^dAf ^ 27rx(Si) = 27r(2 — 2 • genus(St)) ^ An. 


□ 


Using this, Geroch, Jang-Wald discovered a possible approach to prove the Rie- 
mannian Penrose Inequality via the following steps: 

/iSl 

• Let S be the outermost minimal surface in M. Its Hawking mass is y since 
its mean curvature is zero. Notice that this is the lower bound in the Riemannian 
Penrose Inequality. 

• Flow S out by inverse mean curvature flow, and assume that the flow is smooth 
and exists for all time. Let {St} be the flow surfaces. Theorem 2.4 implies that the 
Hawking mass is non-decreasing. 

• Let Sr be the coordinate sphere of radins r in an asymptotically flat coordinate 
chart of M. Theorem 2.3 implies that lim mn^Sr) = mADuiM). 

r^oo 

Here is the npshot: If smooth inverse mean curvature flow of S in M exists for all 
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time, and the flow surfaces approach large coordinate spheres near infinity sufficiently 
fast, then: 


rriADMiM) = lim niHi^t) ^ 

t^OD 



( 2 . 2 . 11 ) 


and this would prove the Riemannian Penrose Inequality. However, smooth solutions 
to inverse mean curvature flow do not always exist. In fact, in the case that the flow 
surface becomes minimal (i.e. mean curvature is zero), the flow is no longer defined 
since the flow speed is the reciprocal of the mean curvature (see (2.2.1)). There are 
other cases where singularities can occur: 


Example 2.2.2 (Inverse Mean Curvature Flow of Disjoint Spheres). Suppose S 
is a disjoint union of two spheres. Inverse mean curvature flow of S will develop 
self-intersection in finite time. 


Example 2.2.3 (Inverse Mean Curvature Flow of Thin Torus). Consider a thin 
torus in obtained as the boundary of an e-neighborhood of a large round circle. 
Thus its mean curvature is positive everywhere. Now starting flowing the torus by 
inverse mean curvature flow (see Figure 2.3). By the first variation of the mean 
curvature (A.f.3) and the parabolic maximum principle, the flow speed has a lower 
bound. As a conseguence the torus will fatten up and eventually the mean curvature 
will become negative in the hole of the torus. Thus, the mean curvature must be zero 
at some point. However, the flow is not defined when the mean curvature is zero. 


2.3 Huisken-Ilmanen’s Approach and Level Set Formulation of In¬ 
verse Mean Curvature Flow with Jumps 


Because inverse mean curvature flow does not always have solutions, Huisken and II- 
manen defined a generalized inverse mean curvature flow which always has solutions. 
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Figure 2.3: Inverse mean curvature flow of a thin torus which develops a singularity 
in hnite time. Picture courtesy of Andrew Goetz. 

The basic idea is that, in this generalized flow, when a surface is enclosed by an¬ 
other surface of less area, it jumps outward to its outermost minimal area enclosure 
(see [25, 5, 6]), and then resume inverse mean curvature flow. Huisken and Ilmanen 
used a level set formulation to characterize this jumping phenomenon. Consider a 
scalar-valued function / on M, and let be the level set of /: 


Et = {xe M\f{x) = t}. 


In this setting, the inverse mean curvature flow equation (2.2.1) becomes: 



(2.3.1) 


In the above, notice that the left hand side is the mean curvature of and the right 
hand side is the reciprocal of the flow speed. Thus when |V/| A 0, equation (2.3.1) 
describes inverse mean curvature flow of the level sets. Note that this formulation 
allows jumps in a natural way, since if / is constant on some region of M, then the 
level sets of / would just jump over that region during the flow. Huisken and Ilmanen 
dehned a weak solution to (2.3.1) using an energy minimization principle, and proved 
existence of such a weak solution by regularizing the degenerate elliptic equation 
(2.3.1). They showed that the Hawking mass is still monotone as in the smooth 
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inverse mean curvature flow case. In particular, the Hawking mass is non-decreasing 
during jumps. In this way, they proved the Riemannian Penrose Inequaltiy (2.1) in 
the case of a single blackhole (i.e. the outermost minimal surface S is connected). 


23 



3 


Inverse Mean Curvature Vector Flow in Spherically 

Symmetric Spacetimes 


In this chapter, we study inverse mean curvature vector flows in spherically sym¬ 
metric spacetimes. The motivation comes from the fact that despite lack of general 
existence theory, inverse mean curvature vector flow always works in spherical sym¬ 
metry. In Section 3.1 and 3.2, terminologies and notations used in later discussions 
are provided. In Section 3.3, we study closed embedded codimension-two surfaces in 
a spacetime and the geometry of their normal bundles. We define mean curvature 
vector fields and Hawking mass. We then define inverse mean curvature vector flow 
of a closed embedded surface. In Section 3.4, we show that spherically symmet¬ 
ric spheres are smooth global solutions to inverse mean curvature vector flow, and 
the Hawking mass is monotonically non-decreasing if the spacetime also satisfies the 
dominant energy condition. 
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3.1 Spacetime, Einstein Equation, Dominant Energy Condition 


A spacetime g, V) considered in this thesis is a connected, smooth, time-oriented, 
four-dimensional manifold with Lorentian metric g of signature (—,-h, -I-,-I-). V is 
the associated Levi-Civita connection. A tangent vector of N is called (see Figure 
3.1) 

1. timelike if g{v,v) < 0; 

2. null if g{v, , n) = 0; 

3. spacelike if g{v, ,n) > 0. 


i 



Figure 3.1: Timelike, null and spacelike vectors in spacetime. 

N is called time-orientable if it admits a smooth timelike vector held T. N is 
time-oriented if such a vector held T is chosen. A vector held X is called future¬ 
pointing if g{X,T) > 0, or past-pointing if g{X,T) < 0. 

(•, •) := (•, ds is used to denote the inner product with respect to g, unless other¬ 
wise specihed. 

A vector held X e F(TiV) is said to be timelike (resp. null or spacelike) if at 
every point p e N, X{p) is timelike (resp. null or spacelike). A submanifold M of N 
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is said to be timelike (resp. null or spaeelike) if every tangent vector of M is timelike 
(resp. null or spacelike). 

Let Ric^ and be the Ricci curvature and scalar curvature of the spacetime 
respectively. The Einstein curvature tensor G is dehned as: 

G := Ric^ - ^R^ (3.1.1) 

We assume that the spacetime satishes the Einstein Equation: 

G = SttT, (3.1.2) 

where T is the stress energy tensor. For any tangent vectors of iV, T{u,v) has 
the physical meaning as the energy density going in the direction of u as observed 
by someone going in the direction of v. The dominant energy condition is: 

T(m, t;)^0, V M, n future-pointing, timelike. (3.1.3) 

3.2 Spacelike Hypersurfaces 


Given a spacetime. Consider a spacelike hypersurface M with globally dehned, 
future-timelike unit normal vector held n on M, and induced Riemannian metric 
g = by restricting the spacetime metric onto M. M is also called a slice in a 
spacetime. Let k be the second fundamental form of M, then the triple {M^,g, k) is 
called a Cauchy data. Given a Cauchy data and the unit normal vector held n, we 
dehne the energy density of M as /i := T(n, n). We can compute that 




1 

Stt 


G(n, n) 


1 

Stt 


^Ric^(n, n) 


^R^ W(n,n)^ 


(By the Einstein equation) 


= ^(Ricqn.„) + iR'') 

= + (trace^A;)^ - \\k\\l) . 


(n is unit time like) 


(3.2.1) 
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where is the scalar curvature of {M,g). We dehne the momentum density of M 
to be a one-form J(-) on M, such that J{X) := T(n, X), for any tangent vector field 
X on M. Then: 

J{-) = T{n, •) = ^ (^Ric^(n, •) - 

= —Ric'^(n, •) (n is normal to M) 

8tt 

1 

= —diVg((A; — trace„A;) ■ g). (3.2.2) 

OTT 

Equation (3.2.1) and (3.2.2) follow from the Gauss and Codazzi equations respec¬ 
tively, and they are called the constraint equations. The dominant energy condition 
(3.1.3) implies that 

^S||J||,. (3.2.3) 

In the time-symmetric case where the second fundamental form k = 0, we see that 
fi = and J = 0. The dominant energy condition thus reduces to R ^ 0. 

3.3 Geometry of Codimension Two Surfaces in Spacetime 

Let S be an closed, embedded, spacelike surface in N with codimension two. We 
assume that S is an oriented surface such that at each point the notion of “outward” 
and “inward” is well-defined. Let g^. be the induced metric on S from the spacetime 
metric g. 

3.3.1 Rank-two Normal Bundle Geometry: Normal Connection and Connection 
One-form 

Let iVS be the rank-two normal bundle of S. Dehne a connection on iVS, denoted as 
V'*", to be the projection of V onto iVS. Notice that since S is spacelike, iVS has an 
induced metric with signature (—, -h). Therefore, at each point p e S, iVpS has four 
quadrants: future-timelike, past-timelike, outward-spacelike and inward-spacelike. 
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Given a local ortlionormal frame { 61 , 62 } of NT,. Suppose 61 is outward-spacelike 
and 62 is future-time like, then the geometry of NpT is depicted by Figure 3.2: 



Figure 3 . 2 : Rank two normal bundle at a point p on a surface with orthonormal 
basis 61 (p) and 62 (p). 

Dehne a linear isomorphism, denoted as “_L”, on each hber NpT of the normal 
bundle as follows: for any orthonormal basis {u, n} of NpT such that u is outward- 
spacelike and v is future-timelike, set := v, and then extend linearly to the entire 
hber. This dehnition is independent of the choice of basis, and is an involutive 
isomorphism. Notice that 

= (v, v} = —1 = —(m, u). 

This isomorphism can be viewed as analogy of the 90° rotation in Euclidean space 
(see [9]). 

Given any tangent vector held X of S. If n is an outward-spacelike unit normal 
vector held along T, then: 

0 = Vx(v,v) = 2(VxV,v). 

Therefore is perpendicular to v, and hence is parallel to v^. From this one can 
dehne a one-form on S, uniquely depends on n, such that: 

a^(X) := <V^n,n^>, VXeF(TS). (3.3.1) 
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This definition yields the following straightforward corollary (see [9]): 

Corollary 3.1. For a smooth section v e r(iVS), the associated connection one-form 
q;„ vanishes if and only if v is parallel with respect to the normal connection V-*-, i.e, 

Vj^v = 0, VXer(TS). 

3.3.2 Mean Curvature Vector Field, Flawking Mass and Inverse Mean Curvature 
Vector Flow 

We define 

If ; TS X TS ^ iVS, (X, Y) ^ (VxY) Itvs (3.3.2) 

to be the second fundamental form of S, where X, Y are tangent vector fields along 
S, and (Vxh^)lvE is the projection of VxV onto the normal bundle of S. Define 
the mean curvature vector field of S to be the trace of the second fundamental form 
with respect to gs: 


Hy. := trace^j,!!. (3.3.3) 

Therefore, Hy is an inward-pointing normal vector field along S. Given any 
normal vector field n, define the mean curvature scalar in the direction of n to be: 

Hn:=-{HY,n). (3.3.4) 

Given (S,5fs) and the mean curvature vector Hy) the Hawking mass of S is 
defined to be: 

mHiV.) := J g{HY,HY)dAY^ , (3.3.5) 

where |S| is the area of the surface S. 
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Definition 3.3.1. Given a spacetime {N,g) and a surface Define the in¬ 

verse mean curvature vector field of S to be: 


Jv := — - 


Hy 




(3.3.6) 


According to our sign convention for Hy:, /s thus defined is outward-pointing. 

Definition 3.3.2 (Smooth Inverse Mean Curvature Vector Flow). Given a closed 
embedded surface Tfi in a spacetime {N^,g). A smooth inverse mean curvature vector 
flow of Yi is a smooth family of surfaces F : S x [0,T] —> N of Y satisfying the 
following evolution eguation: 


■^F{x,s) = TyXx,s), s e [0,T] and (a;,s) e := F(S,s). (3.3.7) 

os 

T > 0 could also be oo. By the first variation of area formula (see Equation A.3.7 
in appendix A.3), the rate of change of area form of the flow surfaces under smooth 
inverse mean curvature vector flow is given by: 


—dAs, = -(Hy,, I-sfidAY^ = dAY,. (3.3.8) 

That is, the rate of the area form of each surface is the area form itself, everywhere 
on each surface. This is a special case of a uniformly area expanding flow first defined 
by H. Bray, J. Jauregui and M. Mars in [10]. 

3.4 Model Spacetime: Spherically Symmetric Spacetime 

In this section, we study spherically symmetric spacetimes. The main motivation 
comes from the fact that, even though inverse mean curvature vector flow lacks a 
general existence theory, smooth solutions still exist in many spherically symmetric 
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spacetimes. Moreover, the Hawking mass is monotonically non-decreasing if the 
spacetime satisfies the dominant energy condition. Thus, it is critical to understand 
the geometry of spherical symmetry. 

Definition 3.4.1. A spacetime {N^,g) is said to be spherically symmetric if its 
isometry group Isom{N'^) contains a subgroup G that is isomorphic to the rotation 
group SO{3); moreover for any point p s N, the orbit of p under the action of G is 
a two-sphere with metric a multiple of the standard round metric. 

From the above dehnition, and standard sphere S'^ share SO(3) as a subgroup 
in their isometry groups, thus share some of the symmetries with S^, and hence 
the term spherically symmetric. 


Proposition 3.2. If is a spherically symmetric spacetime that admits a 

coordinate chart {t, r, 6, 0}, such that g takes the form: 


9 = 


t r 9 

—v‘^{t,r) 0 0 

0 u‘^{t,r) 0 


0 
0 

0 0 0 


0 


0 


(3.4.1) 


0 sin^ 6 j 


where u and v are smooth funetions of t and r only. Then within eaeh t = constant 
slice, smooth inverse mean curvature vector flow of coordinate sphere St^r exists for 
all time with monotonically non-decreasing Hawking mass. 


Remark. Roughly speaking, all spherically symmetric spacetimes outside blackholes 
admit such metrics as in (3.4.1). 

Proof. Let gt^r be the metric on St,r, then 

gt,r = r‘^d9‘^ -\- r^ sin^ 9d(jfl, 
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and its inverse is given by: 


’ sin tl 


Let Ht,r be the mean curvature vector of St^r- Notice that forms a frame 

for the normal bundle of St^r, and thus Ht^ can be computed as follows: 


Ht,r = gU^eAY 


— p, I (Adi dj , dr} ~ 

■*'W <AS,) ‘ (drA) ’ 


m ({^doSeAt) p {VeedeAr)p\ ^ ((^dAAt} p (Vdedp,dr}p 


.. <Va (Ad^AAr). 


<Ad;> + -IaaT " 


- j"(rL% + r«cl) + 0 + + r;,c\) 


is diagonal) 


\ u^J ^ sin^ 0 V ^ 


(See Section A.l) 


-_21^ 

o 

r 


(3.4.2) 


That is 


2 1 

Hf r ~3ir. 

rr rtiL 


(3.4.3) 


The inverse mean curvature vector is: 


^t,r 


Hi A 


-\dr -\3r r 

_ r ' _ _ r ' _ _p\ 

/_2J_^ _2J_^ \ “ AA.,2 “ 2 


(3.4.4) 


Therefore, inverse mean curvature vector flow of St^r is a reparametrization of 


radial flow, and hence is smooth and exits for all time. 
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To prove monotonicity of Hawking mass, recall that: 


mniSt^r) = 


1^ 

IGtt 


1 - 


IGtt 


dAf r , 


(3.4.5) 


St,r 


where dAt^r is the area form of St^r, and \St,r\ is the area of St^r- Note that 

p r* 27 r rn 

= dAt^r = sin 9 dOdcj) 

JSt,r ’ ao Jo 


n 27 r 

Jo 


= 2 (i0 


= 47rr 


Plug the mean curvature vector of St^r (3.4.3) into the Hawking mass equation (3.4.5), 
we get: 

I -^ / 

1 


/47rr^ 

mniS^r) - \l— I 1 


_ 2 1 2 1 
9 ( n^r 1 dAf^>p 


IStt V r r 


til 


1 <■ 4 1 

W^^fi'(dr, dr) dAt^r 


- ( 1 

2 V IGtt 


0 Jo 


2 \ IGtt J^^ ^ u 

^•277 CT7 


4 X 

— — r^sind dO dcj) 


r / 11^ 

2 "" 


r I 1 
-11-^ 


(3.4.6) 


By Equation (3.4.4), the variation of the Hawking mass of St^r along inverse mean 
curvature vector flow is given by: 




V 1-^ + r — 
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Let G be the Einstein cnrvatnre tensor of (iV^, g). G{dt, dt) can be compnted as (see 
Eqnation (A.1.8) in Appendix A.l): 


G{dt,dt) 


2 n ^ o 1 
— 

r U'^ 





(3.4.7) 


By the dominant energy condition (3.1.3), G{dt,dt) ^ 0. Therefore 


1 

2 





+ rGr ^ 0 , 




(3.4.8) 


which implies that It,r{jnH{,St,r)) ^ 0, as desired. □ 

For any spherically symmetric spacetime in Proposition 3.2, it can then be foli¬ 
ated by f = constant spacelike hyperplanes, and each hyperplane can be foliated by 
smooth inverse mean cnrvatnre vector flow spheres (see Fignre 3.3). 



Figure 3.3: Inverse mean cnrvatnre vector flow of coordinate spheres St^r in spher¬ 
ically symmetric spacetime (3.4.1). 
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4 


Spacetimes that Admit Inverse Mean Curvature 

Vector Flow Solutions 


In this chapter we construct non-spherically symmetric, non-static spacetimes that 
admit smooth global solutions to inverse mean curvature vector flows. Surfaces we 
study here satisfy the topological and geometric setups defined in Section 3.3. In 
Section 4.1, we recall the definition of inverse mean curvature vector flow and the 
bad existence theory of such flows. Motivated by the spherically symmetry case, in 
Section 4.2, we construct spacetimes that admit a special coordinate chart (called an 
inverse mean curvature vector flow coordinate chart) in which smooth inverse mean 
curvature vector flow of coordinate spheres exists for all time. We show that there in 
fact exist infinitely many spacetimes that admit such coordinate charts, and hence 
admit smooth solutions to inverse mean curvature vector flow. 

In Section 4.3, we give a coordinate-free analogue of our construction, and show 
that we can actually “steer” a spacetime metric in a certain direction to make the 
mean curvature vector of a surface embedded in a spacelike hypersurface M tangen¬ 
tial to M. Finally in Section 4.4, we discuss some generalizations of the technique 
we use in constructing inverse mean curvature vector flow coordinates. 
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4.1 Motivations from Spherically Symmetric Spacetimes and Main 
Results 

Given a spacetime {N‘^, g) and a closed codimension two surface S with induced 
metric qy. and mean curvature vector Hy, recall from Section 3.3.2 that a smooth 
inverse mean curvature vector flow of S is a normal variation 

F : S X [0, T] ^ M (S, s) - F(S, s) =: 

such that 

^^F{x,s) = TyXx,s), V(a;,s)eS„ (4.1.1) 

where Iy^ is the inverse mean curvature vector of dehned as 


The inverse mean curvature vector flow equation 4.1.1 (same as Equation 3.3.7) is 
a forward-backward parabolic PDE, forward-parabolic in spacelike directions and 
backward-parabolic in timelike directions. Such a PDE lacks a general existence 
theory. However, such PDEs can still have solutions if we start with the “right” 
initial conditions. In spherically symmetric spacetimes, the “right” initial surfaces 
are spherically symmetric spheres. Had we chosen some other sphere to start with, 
inverse mean curvature vector flow is very likely to not exist. This is due to the 
mean curvature vector computation in Equation 3.4.3 in spherical symmetry: it is 
radial and has no components in the timelike direction, therefore the inverse mean 
curvature vector flow of spheres will be contained inside t = constant spacelike slices. 
Geometrically, spherical symmetry eliminates all the timelike “wiggles” of the flow 
surfaces, hence restricting the flow direction to be spacelike. Since the inverse mean 
curvature vector flow equation is backward-parabolic only in timelike directions, in- 
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verse mean curvature vector flow exists without running into singularities in spherical 
symmetry (see Figure 3.3). 

How do we generalize the spherically symmetric case to non-symmetric space- 
times? Note that the mean curvature vector of S is a section of the normal 
bundle iVS, and thus has a timelike component and a spacelike component. Mo¬ 
tivated by the spherically symmetric case, intuitively if the flow surfaces all have 
“purely spacelike” mean curvature vectors, we might hope to get a better existence 
theory. 

Proposition 4.1 ([8]). The family of closed embedded spacelike surfaces is a 
solution to the smooth inverse mean curvature vector flow with spacelike inverse mean 
curvature vector everywhere on the surfaces if and only if there exists a spacelike 
hypersurface a N, such that the mean curvature vector is tangential to M 
at all {x, s), and {Ss} is a solution to the smooth inverse mean curvature flow in M. 

Proof. Given a solution to the smooth inverse mean curvature vector flow with 
spacelike inverse mean curvature vector Js^. Consider the hypersurface M oi N 
dehned by the union of all the surfaces S^, i.e. the “sweep out” region of the the 
flow surfaces. Since the flow is smooth and spacelike, M is a smooth manifold which 
is spacelike as well. ilTs* is tangential to M by the construction of M. Since is of 
codimension one in M, is parallel to the unit outward normal vector held 
along i.e., 

= -Az/s,, (4.1.2) 

at each {x, s) e S* for some smooth positive function A. A is chosen to positive since 
points inward. Therefore 

J _ _ (A I 3') 

(-Az/s,,-Az/s,> A2(z/s,, z/s,> A ffs/ 
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where is the mean curvature scalar of in the direction of z/s^, dehned by 
Equation 3.3.4: 

■= z/E^> = A. 

From equation (4.1.3), we see that indeed is a solution to the smooth inverse 
mean curvature flow in M. 

Conversely, suppose M is a spacelike hypersurface and is a solution to the 
smooth inverse mean curvature flow in M. Assuming that H-£^ is tangential to M 
at each {x,s), we know that is spacelike as well. Moreover, by reversing the 
computations in equation (4.1.3), we have 


Thus is a solution to the smooth inverse mean curvature vector flow equation 
in N, with spacelike inverse mean curvature vectors. □ 

Therefore spacelike smooth inverse mean curvature vector flow solutions in N 
correspond to smooth inverse mean curvature flow solutions in a spacelike hyper¬ 
surface of N with tangential mean curvature vector helds. Huisken and Ilmanen 
dehned a weak notion of inverse mean curvature how in which jumps are allowed. 
This suggests the following dehnition of a weak solution of inverse mean curvature 
vector how: 

Definition 4.1.1 (Weak Solution to Inverse Mean Curvature Vector Flow, [8]). A 
family of spacelike surfaces is said to be a weak solution to the inverse mean 
curvature vector flow equation if there exists a spacelike hypersurface M in N con¬ 
taining Ss such that is tangential to M everywhere for all s e [0, T], and {Ss} is 
a solution to Huisken-Ilmanen inverse mean curvature flow in M (i.e. with jumps). 
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Now we focus on the problem of finding spacetimes with inverse mean curvature 
vector flow solutions. 

Suppose a spacetime {N‘^, g) admits the following special foliation: N is foliated 
by spacelike hyperplanes, and then within each hyperplane, smooth inverse mean cur¬ 
vature vector flow of spheres exists and foliates the entire hyperplane. Consequently, 
the mean curvature vector of the flow spheres are tangential to the hyperplane. If 
N admits such a special “double” foliation (e.g. spherically symmetric spacetimes), 
then N has to be topologically equivalent to x M, where Bi is the closed 

unit ball in 

Suppose admits such a special foliation. We can use this to dehne co¬ 

ordinates that generalize the spherically symmetric coordinates. We define the t- 
coordinate by setting each hyperplane as t = constant, thus the t-coordinate tells us 
which hyperplane we are on. For each inverse mean curvature vector flow sphere, 
dehne A = dvrr^, where A is the area of that sphere. This dehnes a very natu¬ 
ral r-coordinate. Since inverse mean curvature vector how is area expanding, the 
r-coordinate is well-dehned. For simplicity we assume that r ^ 1, i.e. the initial 
spheres on each hyperplane have area 47r. Then dehne {6, 0)-coordinates on an ini¬ 
tial sphere, 0 < 6 < 7i and 0 < 0 < 27r, such that the area form satishes 

dAo = j ^ sm.9ddd(l) = sin 6d6d(l), (4-1.4) 

47r 

where A(0) is the area of the initial sphere. Extend {6, 0) by setting them to be con¬ 
stant along perpendicular directions of the initial sphere, such that {6, 0) coordinates 
are dehned for each sphere. By the extension, ^ will be perpendicular to each sphere. 
The equation for the area form (4.1.4) will be preserved: dAr = sindddd(/)*. See 

* In smooth inverse mean curvature vector flow, -^{dAs) = —{Ig, Hs^dAg = dAg, where dAg is 
the area form of The solution to this equation is dAg = e^dAo. The area of is thus given by 
d.(s) = A{0)e^ = 47rr^, by the definition of the r-coordinate. Thus = A. Therefore the area 

form in the r parameter is dAr = e^dAo = sin 6d9d(j) = sin 6d6d(j) = sin 9d6d(f>. 
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Figure 4.1 for an illustration. 


{N\ 


special 

foliation 



t = 2 


t = 1 


t = 0 


Figure 4.1: Special foliation of a spacetime first foliated by hyperplanes, 

then each hyperplane is foliated by inverse mean curvature vector flow of spheres. 
This generalizes the spherically symmetric case in Figure 3.3. 


Therefore we have proved the “only-if” direction of the following theorem: 


Theorem 4.2. A spacetime {N^,g) is foliated by spacelike hyperplanes, and each 
hyperplane is foliated by smooth inverse mean curvature vector flow of spheres if and 
only if there exists a coordinate chart {t, r, 9, 0} of N, such that in this coordinate 
chart the metric has the form: 


t 



r 9 (j) 

d e f 


9 = 

9 


d 0 0 

e 0 a c 


0 V / 0 c b J 


(4.1.5) 
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where u, v, a, b, c, d, e, / are smooth funetions on N, and the following four conditions 
are satisfied: 


(1) 

o' 

II 

^ 1^ 

(4.1.6) 

(2) 

d d 

dr’ d0 ’ 

(4.1.7) 

(3) 

dAt^r = r^sindddd0 (he. ab — = r"^sin^ d); 

(4.1.8) 

(4) 

Ht^r is tangential to the t = constant hyperplane] 

(4.1.9) 


where dAt^^ and Ht r are the area form and the mean curvature vector of the coordinate 
sphere St,r, respectively. 


Proof of the “if” direction. Given a coordinate chart {t, r, 6*,0} of {N,g) such that 
the g satishes the four conditions, N is then foliated hj t = constant slices which 
are spacelike since the metric has the form (4.1.5). For any t = constant slice, the 
coordinate spheres {St,r} are solutions of a normal flow since (^, = 0. 

We reparametrize the flow by setting s := C* + 2 Inr, where C* is a positive constant. 
Then 


d , , d , ^ ^.ds 2 d 


On the other hand by condition (3) 


(4.1.10) 


— idA) = — (r"^ sin 6dOdcf) = 2r sinddOdcj) = -dA. (4.1.11) 

dr dr r 

Combing the two equations above, we have -^{dA) = dA. Thus, by the hrst variation 
of area formula, {St^r} when reparameterized by = Ge^, are smooth solutions to 
inverse mean curvature flow. By condition (4), the mean curvature vector of St^r 
stays tangential to the slice, therefore {St^r} with the above reparameterization are 
smooth solutions to inverse mean curvature vector flow. □ 
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Definition 4.1.2 (Inverse Mean Curvature Vector Flow Coordinate Chart). If a 
spacetime admits a coordinate chart {t,r, 6^,0} such that the four conditions 

( 4 - 1 - 6 ), ( 4 .1.7), (4-1-8) and (4-I-9) are satisfied, then {t,r, 6*,0} is called an inverse 
mean curvature vector flow coordinate chart, and N is called a spacetime that admits 
an inverse mean curvature vector flow coordinate chart. 

We sometimes refer the fourth condition (4.1.9) as the steering condition, as it 
forces the coordinate spheres to stay inside the spacelike hyperplane during inverse 
mean curvature vector flow. 

Many spherically symmetric spacetimes admit an inverse mean curvature vector 
flow coordinate chart (e.g. coordinate chart (3.4.1) with d = e = f = c = 0, and 
a = r‘^,b = r'^ sin^ 9 and radial mean curvature vector by Equation (3.4.3)). However, 
given an arbitrary spacetime it is generally impossible to reparametrize it 

with an inverse mean curvature vector flow coordinate chart (e.g. a spacetime that 
is not topologically equivalent to (M^\i?i) x M). However, is it possible to construct 
a spacetime that admits such a coordinate chart? In the next section we show that 
we can actually construct many such spacetimes: 

Proposition 4.3 (Existence of Spacetimes That Admit an Inverse Mean Curvature 
Vector Flow Coordinate Chart). Let U := x M c There exist infinitely 

many spacetime metrics of the form of (4.1.5) that admits an inverse mean curvature 
vector flow coordinate chart. 

Combining Proposition 4.3 and Theorem 4.2, we have the following main theorem 
of this thesis: 

Theorem 4.4 (Main Theorem). There exist infinitely many non-spherically sym¬ 
metric, non-static spacetimes that admit inverse mean curvature vector flow coordi¬ 
nate charts. Given such a spacetime U with an inverse mean curvature vector flow 
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coordinate chart (t, r, 9, 0) and the constructed spacetime metric g. The coordinate 
spheres St^r contained in each t = constant slice, when reparameterized by r^ = e^, 
are smooth global solutions to the inverse mean curvature vector flow equation. 

4.2 Construction of Spacetimes That Admit Inverse Mean Curvature 
Vector Flow Coordinate Charts 

In this section we prove Proposition 4.3. Let U = x M. It is easy to construct 

a spacetime metric g that admits a coordinate chart {t, r, 6, 0} that satishes condition 
(4.1.6) and (4.1.7). Simply dehne 


t 



r 9 (f) 

d e f\ 


g ■= 

9 


d 0 0 

e 0 a c 


<t>\ f 0 c b J 


(4.2.1) 


where a, b, c, d, e, /, u, v are arbitrary smooth functions on U. Choosing two of the 
three variables a, b and c such that ab — c^ = r^siM 9 satishes condition (4.1.8). 

The fourth condition requires Ht j., the mean curvature vector held of St^r, to be 
tangential to the t = constant slice. This is equivalent to requiring Ht^r to be parallel 
to We compute the conditions on the metric components such that this is true. 

Lemma 4.5. The determinant of the spacetime metric g in (4-2.1) is given by: 

\g\ := deflg) = (—— df){ab — C) + eu^{cf — be) + fu‘^{ce — af) (4.2.2) 
= (—— d^){ab — C) + u^{2cef — be^ — ap) (4.2.3) 
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Moreover, the eoordinate representation of the inverse {g is given by: 


t 


r 




t 

r 


v?{ab—c^) 

Isl 

—d(ah—c^) 


~-d{ab—c^) 

—v^{ab—c^)-\-f{ce—af)-\-e(cf—be) 

m 


u'^(cf-be) 

\a\ 

-d{cf-be) 

\a\ 


(ce—af) 

\a\ 

—d(ce—af) 


e 


u^jcf-be) 

lal 


-d{cf-be) 

Isl 


—v?v'^b—v?f‘^—bd^ 

Isl 


I u^(ce—af) 


-d{ce—af) 


u^v‘^c-\-u^ef-\-cd^ 


\a\ 


Proof. See Section A.2.1 in Appendix A. 

^.2.1 Geometry of St^r and the Normal Bundle NSt^r 


v?‘v^c-\-v?ef-\-cd^ 

W\ 


—uva—ue 


1^1 


7 

(4.2.4) 

□ 


Fix a. t = constant slice M. Let St^r be a coordinate sphere in M. We endow St^r 
with the induced metric from g, denoted as gs- Then in the {6, 0} coordinate system, 
gs has the following representation: 


9s ■= g\st,r = 


9 (f 

9 I a c 

c b 


(4.2.5) 


Thus its inverse metric is: 


9 s' 


9 (p 
1 9 I b -c 


ab — 


—c a 


(4.2.6) 


The normal bundle NSt^r of St^r is of rank-two. ^ is a nonzero section of NSt^r, 
and thus can be used as a basis for the normal bundle. Let n be a complementary 
basis vector held of the normal bundle that is orthogonal to Since ^ is out¬ 
ward spacelike, we can assume that n is future timelike. Therefore using the basis 
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5 5 5 5 

n = + X— + y— + z—, 

dt dr ^ do d<p 


with X, y, z yet to be determined, such that 


i> = 0; 


(n, i> = 0; 

(n, A) = 0. 


We have three equations and three unknowns which give us 


5 —dd cf — be d ce — af d 

^ 5t 5r a6 — 55 ab — 50 ’ 


(4.2.7) 


Lemma 4.6. 


,n, n) = 


det{g) 


det{g) 


u^{ab — c^) u^det{gs) 


(4.2.8) 


where we set l^'s'l := det{gs)- 

Proof. See Section A.2.2 in Appendix A. 

Remark. Caution that since n is timelike, (n, n) < 0. Thus 


M9 = {-(n,n)Y'^ = 


ITTI \ 1/2 


u'^lOsl 


(4.2.9) 


Let {er,en} be the normalized orthonormal frame obtained from {^,n}: 


5 ,II 5 II 15 n 

c I Is ~ ) Cn •= 7] j7 

dr dr u dr n 


ht,r 
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be the outward-pointing inverse mean curvature vector. The geometry of the normal 
bundle of St^r is given by Figure 4.2 below. Recall that Ht^r points inward by our 
convention. 


fltr 

{Hur. H,.r) 




NSt,r 

Figure 4.2: Normal bundle of coordinate sphere St^r- 

4-2.2 Mean Curvature Vector of St,r 

Now we have a coordinate sphere {St,r,gs) with induced metric gs in (f/, g,'V), where 
V is the Levi-Civita connection with respect to g. We have an orthonormal frame 
{er,en} of the normal bundle NSt,r dehned in the previous subsection. In this sub¬ 
section we compute the mean curvature vector Ht,r of St,r- Recall the dehnition of 
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Ht r from Equation (3.3.3): 

Ht^r = trace^gll 


ott 1 A 




dx'' 


NSt,r 


(By definition of Ht^r 


< O _ O _ p I -fi -fit'l 
dx^ ’ dxi t 50 ) d(j> > ' 


(By definition of II(^, ^)) 




({er,e„} is an orthonormal frame of NSt^r 



ox^ 

(4.2.10) 

We compute the two parts in Ht y separately. First of all 



\ ^ O/V7 ^ ^ \ ^ *i/F^ ^ ^ \ 



1 t 

= -^5 ^ij9tr + ^ij9rr) {t IS perpeudlcular to 


- r‘„ + 29 ''»r‘^ + 

+ + 2fe% + 9yi^«) 

lip - 1 - f - f - Y- 

— .y [9tri^^ee 2crg0 + aB^^) + 

- 2cr;^ + ar;^)] 


” « ISsI *’'’■^**^1 

(4.2.11) 

where 

(.) := bTl, - 2cr‘^ + oF;^, 

(4.2.12) 

and 

(★★) := — 2crg^ + aF^^. 

(4.2.13) 


We now compute all the related Christoffel symbols. 
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^ee ~ 9 {9^^i‘^9et,e 900,t) + 90$,r) + 9^^{‘^900,0 900,0) + 9^^{‘^904>,0 


900, 

1 


~ 2 5'6ie,t) 9^^900,r + 9^^900,0 + 9^’^{‘^90<j,,0 900,<j,)^ 

= ^p| — c^)(2e,0 — a^t) + d{ab — c^)a,r + u^{cf — be)a,0 + u^{ce — af)(2c,0 

(4.2.14) 


— 2 + 90t,<t> ~ 90(j,,t) + 9^^{9<t>r,0 + 90r,<j> ~ 90(j>,r) + 9^^{9<j,0,0 + 900,4> ~ 90(f>,0) 

9 ^9(p(p,0 4 “ 90(p,(p 90(i},(p)'^ 

= \ [f\9<pt,0 + 90t,4> - 904>,t) - t'’90<p,r + f^900,<p + f^94><p,0) 

= 777117 (u‘^{ab — c^){f,0 + e^0 — Cj) + (i(a6 — c^)cr + u^{cf — be)a^^ + M^(ce — af)b^g) 
^\9\ ^ ' 

(4.2.15) 


^ - 94 > 4 >,t) + 9 ^''{‘ 29 ^r,<p - 1/00,r) + l/*^(2^00,0 - 1/00,0) + 9 ^^{‘ 29 ^^,<p 

= 2 {9^\‘^9^t,<f> ~ 9^^,t) ~ 9^''9<p4>,r + ^*^(2^00,0 — S' 00 , 0 ) + ^*'^^ 00 , 0 ) 

= ^ (^M^(a6 - c^)(2/,0 - 6,t) + (i(a& - c^)&,,. + u^{cf - be){ 2 c,^ - 6,0) 

+ u‘^{ce — af)b,^ (4.2.16) 



48 




9ee,t) + 


9ee,r) + 9^\‘^9e0fi “ 90e,e) + 9''^{‘^90^fi 


— 2 {9^^{‘^90t,0 900,t) 9^^900,r + 9^^900,0 + 9^‘^{‘^904,,0 900,(jSj 

= :;^7ir| f — d{ab — (^){2e,0 — a,t) — \—v^{ab — c^) + f{ce — af) + e{cf — he)\a,r 

A9\ ^ 

+ {-d){cf - he)a,0 + {-d){ce - a/)( 2 ce - a,^)'j ( 4 . 2 . 17 ) 


^ 0 ij> 2 (.9(f)t,0 4" 90t,<f> 90(f),t) 4” 9 i.9if>r,0 ”4 90r,if) 90(f),r) 4“ 9 i.9<i>0,0 4“ 900,(f} 90<f>,0^ 

4~ 9 ^9<ff(p,0 4“ 90(i},(i} 90(p,(p)^ 

= \ + 90t,4> - 904>,t) - f''90^,r + 9"'^900,4> + 9^''^9 4 ,^, 0 ) 

= ^ ( - d{ab - c^){f ,0 + e ,0 - c*) - [-w^(a 6 - c^) + /(ce - af) + e{cf - 6 e)]c^ 

+ {—d){cf — be)a^fj) + {—d)(ce — af)b^g^ (4.2.18) 


^ J + 9''''{‘^9<t>r,^ - 1/00,r) 4 - 9^\‘^9^0,4> - 9^^,0) + ^''‘^(2^00,0 

= 2 {^9~ 9<f>4>,t) ~ 9^''94>4>,r + ^^^(2^00,0 “ ^ 00 , 0 ) 4- 

= ^ - c?(a& - c^)(2/,0 - b,t) - [-v^{ab - c^) + f{ce - af) + e(c/ - 6e)]6,^ 

+ {-d){cf - 6e)(2c0 - 6,0) + {-d){ce - af)b,^ (4.2.19) 
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With all the Christoffel symbols computed, we have: 

(.) = (-r‘, - 2crl + arl, 

= - c^){2befi - ^ - 2ce,^ - 2cf,e + 2^ + 2af,^ - 

+ d{ab — c^)[a^rb — 2cc^r + a&,r] + u^{cf — be)[a^gb — 2a,(/,c + 2aC(/, — ab^g] 

+ u^{ce — af)y2bcfi — a^^b — 2b (4.2.20) 

The terms cancel in the second line above since ab — = r* sin^ 9, and thus is not a 

function of t. 

Next we compute 

(**) = bVlg - 2cr0^ + aTl^ 

= - d{ab - c^)[2be^g - ^ - 2ce^^ - 2c/,^ + ^ce/ + 2a 

— [—v^{ab — c^) + /(ce — af) + e{cf — 6e)][a,r& — 2cCr + a&,r] 

— d{cf — be)[a^gb — 2a,+ 2ac^^ — ab^g] — d{ce — af)[2bc^g — a^^b — 2b^gc + a6,<^] 

(4.2.21) 
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Now plug (4.2.20) and (4.2.21) back into (4.2.11): 


'u^{ab — c^) — u^d{ab — c^)][26e^e — 2ce^^ — 2cf^g + 20 /^^] 


_ 1 1 1 

u l^sl 2\g\ 

+ {d^{ab — c^) + u^v^{ab — c^) — fu^{ce — af) — eu^{cf — 6e)][a,r& — 2cc^r + a&,r] 
+ \du^{cf — be) — u^d{cf — be)][a^gb — 2a + 2ac,<j(, — ab^) 

+ \du^{ce — af) — u^d{ce — a/)][26c0 — a^^b — 26 gc + ab^^ 

1 1 1 


u l^sl 2\g\ 
+ ab^r\ 


— \d^{ab — c^) + u^v^{ab — c^) — fu{ce — af) — ew{cf — be))\a^rb — 2cc^ 


1 1 1 


-(-|^|)(a6-c ),, 


u l^sl 2|^| 

1 4r^.s«r^ 1 2 1 


u 


r"^.s«r^ 2 


r u 


(By Equation 4.2.2) 
(4.2.22) 


Next we compute the second term in (4.2.10): 


|n||g-"“ ■ ■"'dt' 


^5 rL<n, n> 


n 


(n — ^ is spacelike) 


= -gfT\^M\g = -WnUds^^ee + + 9t%7 








n 


l^sl 


■w 


(4.2.23) 


where (★) is computed in (4.2.20). Now plug (4.2.8) into the above, we get: 


9|<Vi/-,e„> = 

Ba:^ OX^ 


- —'] ' —(*) = --^^^(*). 

u'^\9s\) Ifi-sr u \gs7^ 


(4.2.24) 


Combing (4.2.22) and (4.2.24), we have: 
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Proposition 4 . 7 . The mean eurvature veetor Ht ^ of coordinate sphere St^r in the 
spacetime metric g (4-2.1) is given by: 

Ht,r = ---Or + (4.2.25) 

ru u \gs\ ' 

where (★) is given by (4-2.20). 

Corollary 4 . 8 . The mean curvature vector Ht r of coordinate sphere St^r is parallel 
to ^ everywhere if and only if {*) = 0, where (★) is given by ( 4 - 2 . 20 ). Moreover, in 
this case, the mean curvature vector eguals: 


^ 2 1 

Hfj. = - Cr- (4.2.26) 

’ r u 

Notice the similarity between the mean curvature vector expression above and 
the mean curvature vector in the spherically symmetric case in (3.4.3). 

Proof. This is quite straightforward since Ht^r is parallel to ^ everywhere if and only 
if On) = 0. Expand this out we get: 


0 (Ht^j.,en) \ 1 ^ (*)On,On/ 




u \gs\^i‘^ 


Since u ^ 0 , the above is equivalent to 


(*) = 0 , 


as desired. □ 

Therefore the fourth condition in the dehnition of inverse mean curvature vector 
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flow coordinate chart is equivalent to (★) = 0, that is: 

0 = “ c^)[ 2 &e ,0 - 2 ce,^ - 2 c/,e + 2 a/, 

+ d{ab — c^){ab — c^),^ + u^{cf — be)\a^eb — + 2 ac,^ — a 6 ,e] 

+ u^{ce — af)\2bcfi — a^^b — 26,^0 + a 6 ,(^]| 

Or equivalently, using ab — sin^ 

4?^^ sin^ 0 

0 = [26e,e - 2 ce,</, - 2c + 20 /,^] + d -^- 

cf — 6c CC — Qjf 

+ . ■ 2n [^,o^ ~ + 2ac0 - ab^g] + . o^ P&cg - a^^b - 2b+ a6,^]. 

sin 6 r* sin tl 

(4.2.27) 

Here is the upshot: Equation (4.2.27) is zeroth order in the metric component 
d, thus we can choose two of the three variables a, b, c, and e, /, u, v all together 6 
variables, and solve for d explicitly: 

( cf — be 

^ ~ , o . 2 n \ - 2 ce,<^ - 2 c/,^ + 2 a/,<^] + . . o ~ 2 a,^c + 2 ac 0 - ab^g] 

4r'^ sin 6 y sin 6 

Qp — n f 

+ , . - a,4>b - 2b^gc + ab^^] [. (4.2.28) 

sm 0 ) 

Thus all four conditions (4.1.6), (4.1.7) and (4.1.8) and (4.1.9) are solvable with 
inhnitely many solutions. Combing the above, we have proved Proposition 4.3 and 
hence the main theorem 4.4. 

There are six degrees of freedom in constructing our spacetime metric g that 
admits an inverse mean curvature vector flow coordinate chart. Moreover, the six 
free variables do not need to be spherically symmetric. The spherically symmetric 
(3.4.1) metric is a special case of this large set of spacetime metrics. It is unknown if 
perturbations of spherically symmetric spacetime with inverse mean curvature vector 
flow coordinate chart still have inverse mean curvature vector flow solutions. More 
specihcally, we conjecture that: 
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Conjecture 4.1. Given Minkowski space with inverse mean curvature vector flow 
coordinate chart that can be smoothly extended to the boundary, consider a pertur¬ 
bation of the spacetime metric. The resulting spacetime still admits inverse mean 
curvature vector flow solutions (in a single spacelike hypersurface) that exist for all 
time. 

Notice that in Equation (4.2.28) if all the variables are smooth then d will be 
smooth except possibly when sin6' = 0, since d is not dehned by our formula there 
(see (4.2.28)). This happens at the north {9 = 0) and south pole {9 = n), which are 
two coordinate chart singularities, not metric singularities of the spacetime. 

If c, e, / are chosen to be 0 at a neighborhood of the north and the south pole, 
then the right hand side of (4.2.28) will be zero there. In this way d can be extended 
smoothly across the two coordinate chart singularities, and will be smooth on the 
entire spacetime. 

Another way to extend d smoothly over the coordinate chart singularities is to 
choose the metric to be spherically symmetric in a neighborhood of the north and 
south pole: 


a = 


t r 9 

-v‘^{t,r) 0 0 

0 u^{t,r) 0 


0 

0 


0 

0 


0 
0 
0 


(4.2.29) 


0 r^ sin 9 j 


This has the advantage that d = 0 around the coordinate chart singularities. Then 
extend a, b, c, e, / smoothly to the entire spacetime while maintaining the condition 
that ab — = r^sirfl 9. The resulting metric still satishes the four conditions since 

smooth inverse mean curvature vector flow of spheres exists for all time in spherically 
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symmetry, d will be smooth since d = 0 identically. 

One could study more general asymptotic conditions for (4.2.28) to be smooth 
and bounded as 9 approaches 0 or n, but we choose not to discuss it further here. 

Remark. We can actually prove that Ht^r takes the form (4-2.26) in inverse mean 
curvature vector flow coordinates without computing it out explicitly. We now show 
that as a sanity check of our computation. In inverse mean curvature vector flow 
coordinates Ht^r is parallel to thus let A = \{t,r,9,(l)) such that 

Ht,r = A^. 

or 

The inverse mean curvature vector is now 


It,r 




t,r 




1 d 




By the first variation of area formula (A.3.7), the rate of change of the area form of 
St^r under outward radial flow is given by 

fdAs,, = -</?,„ f)dAs,„ 

d 


= —(Htr, AnOdOdf) 
’ dr 

= —Xu'^r'^ sin 6dOdcj) 


(4.2.30) 


Notice the left hand side of the above eguals to 


—dAs, = — (r^ sin dddd^) = 2r sin 
dr ’ dr 


Thus matching the two sides we get: 


2r sin9d9d(j) = — Xmr^ sin 9d9d(j), 
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that is \ Therefore 

r J 




2 l d _ 2 1 

r dr ru 


which is the same as (4-2.25). 

4.3 Coordinate Free Analogue and Steering Parameters 


In the previous section we have shown that there exist many spacetimes that ad¬ 
mit inverse mean curvature vector flow coordinate chart, in which the coordinates 
spheres are solutions to the inverse mean curvature vector flow equation. The fourth 
condition (4.1.9) in the definition of inverse mean curvature vector flow coordinates 
can be viewed as a steering condition that keeps the flow direction of coordinate 
spheres tangential to a spacelike hypersurface. 

Given a spacetime V), a spacelike hypersurface (M^, g) with induced met¬ 

ric g, and a closed embedded surface (S, g-^) in M with induced metric gs- Assuming 
the normal bundle of S is trivial, there is a unique unit outward normal vector held 
of S in M, denoted as e^. Let be the unit outward normal vector held of S in 
that is perpendicular to e^. Since M is spacelike, Cr is spacelike and e* is timelike. 
Define a local coordinate chart on S, and let {^, be the local coordinate 

frame. Figure 4.3 depicts the above set up, in which the mean curvature vector Hs is 
not necessarily tangential to M. We can extend the frame {e^, e^, to a frame 

on a neighborhood of S in N, and we will identify the frame with its extension. 

Since the local frame { 64 , 6 ^, is not fully a coordinate frame, the commu¬ 

tator coefficients O’), defined as: 

d d 

[ttj, Oji] C^jOtk-i oii G {e^, Cf, } (4.3.1) 
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are not necessarily zero. We need these coefficients to compute the connection coef- 
hcients with respect this frame later. 





Figure 4.3: Setup of inverse mean curvature vector flow steering. 


With the local frame {e*, e^, of the tangent bundle of N, let the associated 

dual frame of the cotangent bundle of N be {/3t, /3r,d6,d(j)}, where 

= 6ij, ijE{t,r}. 

With respect to this dual frame, we can write the spacetime metric g as 


et 



0 


e_ 

ee 

0 


m 


d_ 

89 

dip 


0 

V 0 


1 0 

0 a 

0 c 


8 _ 

dip 



0 


b J 


(4.3.2) 


where a, 6, c are smooth functions on N. 

We want to change the metric so that H^, is tangential to M. Recall that the 
fourth condition (4.1.9) in the construction of the inverse mean curvature vector 
flow coordinates is a zeroth order equation for d, the (t, r)-metric component. This 
motivates the following dehnition: 
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Definition 4.3.1 (Steering of Spacetime Metric). Given a spacetime metric g. A 
metric Qq on N is called a steering ofg if 

gQ-=g + Q{Pt®Pr + /dr®Pt) (4.3.3) 

for some smooth function Q e C°°(iV). Q is called a steering parameter. 

Note that the coefficient matrix for the steered metric gQ is given by: 


et 



Q 


e_ 

ee 

0 


8_ 

d4> 



[gq] = 


e_ 

ee 


Q 

0 


0 


1 0 

0 a 

0 c 


0 


bj 


(4.3.4) 


Thus geometrically, et and are not necessarily orthogonal to each other in the 
steered metric (see Figure 4.4). 



Figure 4.4: After steering the spacetime metric g, et and are not necessarily 
orthogonal to each other. The metric g on M remains the same. 


Definition 4.3.2 (Area Expanding Condition). S c M given as above is said to be 
area expanding if 


er{ab — (?) > 0 . 


( 4 . 3 . 5 ) 







Theorem 4.9. Let a spacelike hypersurface and a closed embedded sur¬ 

face (S^, 5 fs) in M be given as above. If Yi is area expanding, then there exists a 
unigue smooth steering parameter Q = Qs ^ C°°(iV), such that in the steered space- 
time metric gq, is tangential to M everywhere on S. 


Proof. To show that there exists Q such that H^, is tangential to M, it suffices to 
hnd Q such that H^, (with respect to gq) is parallel to e^. 

Notation: In the following the subscript Q will be dropped for simplicity and all 
the spacetime metric will be referring to the steered metric gq. 

The computations are similar to the inverse mean curvature vector flow coordi¬ 
nates case. We divide the computations into hve steps: 

Step 1: Pick a local normal variation of S and let be the variational surfaces, 
s e (0, e). Since Sq = S is area expanding, we can assume that this normal variation 
is also area expanding. Extend {et, Cr} to a neighborhood of S in iV such that 
remains outward unit normal to each S* in M. Dehne an r-coordinate by requiring 
the area of Yg to be ^(<s) =: dvrr^. Since the variation is area expanding, r is well- 
dehned. Extend to a neighborhood of S in iV such that they are constant 

along normal directions of S in M. By this extension ^ is perpendicular to each Yg. 
Thus there exists a function A such that = A^. 


d d d d 


(4.3.6) 


This implies that C% = 0. Similarly = 0. 

Step 2: After steering the metric g, and et are not necessarily orthogonal to 
each other. Let n be the normal vector of S such that (n, ef) = 0. The computation 
of n is the same as the inverse mean curvature vector flow coordinate case (4.2.7), 
and we obtain: 

n = et — Qcr. (4.3.7) 
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Note that (n, n) = (et — Qe^, et — Qe^') = —1 — 2(5^ + = — (1 + Q^), which agrees 

with (4.2.8) with our metric Qq. 

Let Cn ■= be the unit timelike normal vector. 

Step 3: With respect to the orthonormal frame {e„, e^} of the normal bundle iVS, 
we can write the mean curvature vector as 

Ht, = {Hs, er)er — {Hs, en)en- (4.3.8) 

Therefore Hs is parallel to Cr if and only if (iLs, Cn) = 0. Now we compute the 
condition on Q such that this inner product vanishes. 

(Hj:, en) = (traceg^n, e„> = c/g<(VaTj) |jvs, e^) (d^, dj e ^}) 

= g'l {(ydidj, er}er - <VaTj, e„>e„, e„) 

~ 9T,dySidji(ird)(Gm(ird) = di^j ■, Gri) (^Cn; = 1) 

= j^^aiydidj^n) = n> 

= i^4H(-(i + «^)) - 

= -(1 + QY^a'H ( 4 . 3 . 9 ) 

where u is the connection coefficients of the connection Vq with respect to the 
metric gg and the local frame {et,er, which is not necessarily a coordinate 

frame. This is an important difference from the inverse mean curvature vector flow 
coordinates case. 

Recall that 

+ 9ji,i - 9ij,i + Ciij + Ciji - Ciji), (4.3.10) 

where Ciji := giy^, and C™ is the commutator coefficients dehned in (4.3.1). 
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Notation : we use to denote ek{gij), not necessarily a coordinate derivative. 
Step 4: Now we compute the four connection coefficients as follows: 

^ee = ~ 9ee,i + - Ceei) 

1 _. _ _ _ . _ 

= ^9^\‘^9ei,e - 9ee,i + '^9ejCle - is a coordinate frame) 

= 2 “ ~9de,t + ‘^9ejCie) + - 9ee,r + ‘^9ej99ig) + 0 + o| 

= 2{9^^^~9e9,t + ‘^9eeCw + ‘^9e<j,Cte) + 9^''{-9ee,r + ‘^9eeC% + ‘29d4>Cte)"\ 

= 2 |^|—+ 20(7)^) — Q{—a^r + 20(7^^ + 20(7^0)1 (4.3.11) 




1 

2 S' + 94>i,e ~ 9<i,e,i + C'itAe* + C'je-/- “ C'fliei) 


-^ 9 ^^ 901 , 4 ) + + fi'ejC'i!/. + “ & 



(4.3.13) 
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0 


2^ ^9(j)i,(j) 9(f)(j),i ‘^9(j>j^i(p ) 

= 2 ~ ‘^94)0^14, + ‘^94>4>Cf^) + “ ^00,r + ‘^ 94 , 0 ^% 

+ 2^00C^0)| 

= + 2<0 + ‘^bCf^) - Qi-br + 2cC% + 2bCt^)] (4.3.14) 

Step 5: Now plug (4.3.11), (4.3.12), (4.3.13) and (4.3.14) back into (4.3.9), we 

get: 

(Hj:, en) = -(1 + \ ^2 (b^o0 - ci^0<p “ cwj-e + 0 ^ 00 ) 

= -(1 + \ ^ 2 ^bulg - 2ajl^ + au^^) (By (4.3.13)) 

= —(1 + Q‘^y ^'^^2 2\g\ ‘^o,bC^Q + 2})eCfg + 2cc^t — ‘^c^C% — 

- 2cteq^ - 2^Ct^ - ab^t + 2cieC^ + 2abCf^) - Q{-a,rb + 2abC% + 

+ 2ccr - ‘2c^C% - - 2ax<^ - 2c^Cf^ - ab^r + 2cteC:^ + ^abCf^) | 

= —I - et{ab - c^) + 2{ab - c^){C^g + Cf^) 

-Q (^-e^(a6-c2)+ 2(a6-c2)jC:^^Ke^ 

= - ^ —|e,.(a6 - c^)Q - et{ab - (?) + 2(a6 - ?){C% + Cf^)^ (4.3.15) 

since C% = Cf^ = 0 by Equation (4.3.6). Therefore Hy. is parallel to if and only if 
er{ab — ?)Q — et{ab — ?) + 2{ab — ?){C^g + Cf^) = 0. (4.3.16) 
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Notice that Equation (4.3.16) is zeroth order in Q. Since S is area expanding, 
er{ab — c^) > 0, and hence we get a unique solution 


Q 


etjab - c^) - 2{ab - c^)(C'fg + Cf^) 
er{ab — c^) 


(4.3.17) 

□ 


Lemma 4.10. e^iab — c^) = 0 if and only if = 0, i.e. T, is a 

minimal surface in M, where He^ is the mean curvature scalar of S in the direction 
of Cr. 


Proof. For the first claim, note that 

er{dAYf) = er{\lab - c^jdefd^) = ^ ^ ^ ^ er(ab-c^)/30/d^ = ^^y^j^^—^er(ab-c^)dAs. 


2\fah 


2(ab — c^) 


(4.3.18) 


On the other hand by the hrst variation of area formula (A.3.7), 


er{dAY,) = er)dAY, = He^dA-£. 


(4.3.19) 


Combing the two equations we get: 

2i7er ■ {oh — c^) = er{ab — c^). (4.3.20) 

Therefore er{ab — c^) = 0 if and only if He^ =0. □ 

An application of Theorem 4.9 is to generate more examples of inverse mean 
curvature vector flow solutions. Let (A^"^,^) be a spacetime and {M,g) a spacelike 
hypersnrface with induced metric g. Suppose is a solution to the smooth inverse 
mean curvature flow in M. Let Hg and dAg be the mean curvatnre vector held and 
area form of S^, respectively. By the hrst variation formula, 

^{dAg) = dAg. (4.3.21) 
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Theorem 4.9 allows us to steer the metric smoothly along inverse mean curvature 
flow to keep Hg tangential to M. Note that (4.3.21) still holds after the steering since 
we are not changing the metric on M. Therefore by Proposition 4.1, in the steered 
metric {Ss} is a solution to smooth inverse mean curvature vector flow equation. 

One could generalize this technique to study weak solutions to inverse mean 
curvature vector flow equation (defined in (4.1.1)) using solutions to Huisken-Ilmanen 
inverse mean curvature flow (with jumps) in a hypersurface, but we will not give a 
rigorous treatment here. 

4.4 Generalizations 


Given a surface {Ti,g-s) inside a spacetime {N^,g). Let I be the inverse mean curva¬ 
ture vector. A more general flow than inverse mean curvature vector flow is to flow 
out S in the following direction: 

+ (4.4.1) 


where the T operation is a linear isomorphism on the normal bundle defined in 
Section 3.3.1, and /d is a constant on each flow surface (hence is only a function of 
the flow parameter) such that — 1 < /d < 1. 

Therefore inverse mean curvature vector flow corresponds to the case where /3 = 0. 
The procedure for constructing spacetimes with inverse mean curvature vector flow 
solutions can be generalized to constructing spacetimes in which this more general 
flow exists. The idea is to construct a spacetime metric g that admit a coordinate 
chart {t, r, 9, 0}, such that conditions (4.1.6), (4.1.7), (4.1.8) and the fourth condition: 


^ is parallel to —. 

dr 


(4.4.2) 


are satisfied. 
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5 


Uniformly Area Expanding Straight Out Flows 


In this chapter we study uniformly area expanding straight out flows, or simply 
straight out flows of spacelike surfaces in a spacetime. Consider a spacetime (W,^) 
and a closed embedded surfaces S with the induced metric gj]. A straight out direc¬ 
tion of S, hrst studied by M. Mars, E. Malec and Simon [30], is a normal vector held 
that has “minimal variations” along S. Such a normal vector held is a minimizer 
of a natural functional dehned on the normal bundle. The Hawking mass is also 
monotonically non-decreasing under smooth straight out hows assuming the space- 
time satishes the dominant energy condition. A condition for a spacetime to admit 
straight out how coordinate charts is derived in this chapter. Complete understand¬ 
ing of such spacetimes is still work in progress. 

5.1 Background and Straight Out Flow Coordinate Chart 

Let be a time-oriented spacetime and (Ti‘^,g-£) be a closed embedded sur¬ 

face. The normal bundle of E, XS, has an induced metric of signature (—, -I-). On 
each hber NpT,, the nonzero vectors get partitioned into four quadrants: outward- 
spacelike, inward-spacelike, future-timelike and past-timelike. Let denote 
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the subbundle of iVS that consists of outward-spacelike normal vector fields of unit 
length. Given a smooth section v of its associated connection one-form ai^ 

on S is defined by: 

a,(X) := <Viu,u^> 3 , VX e r(TS). (5.1.1) 

where _L is a hberwise linear isomorphism dehned in Section 3.3. Notice that {z/, 
forms an orthonormal frame of iVS. 

Given another smooth section V of there exists a constant 0 > 0 such 

that 

V = cosh 6v -I- sinh 

= smhOv -I- (5.1.2) 

Geometrically {z7, p-*-} can be viewed as hyperbolic rotation of {z/, by angle 9. 
The associated connection one-forms are related by: 

Lemma 5.1 ([9]). Let he the assoeiated eonnection one-form ofV, then ajj is 
related to a^, by 

a-jj = — dO. (5.1.3) 

Let E be an energy functional on U^N(T,) such that: 

E(z/) := J llV^z/fdAs, Vz/er(G+iV(S)). (5.1.4) 

H. Bray and J. Jauregui proved the following proposition: 

Proposition 5.2 ([9]). z/ e r(iVE), outward spacelike and is of unit length, is a 
minimizer of E if and only if diusi^au) = 0. 

Intuitively, a minimizer of i? is a normal vector held with “minimal variations” 
along E. Bray and Jauregui also proved the existence of such minimizers: 
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Proposition 5.3 ([9]). Minimizers of E exist, and any two such minimizers differ 
by a hyperbolic rotation as in (5.1.2) of a constant angle 9. 

Proof. We present the proof in [9]. Fix a section u of U'^N(T,). Given any other 
section z7, z7 and u are related by a hyperbolic rotation as in (5.1.2). Let a,y and Ou 
be the associated connection one-forms of u and V, respectively. By Lemma 5.1, aj 7 
is a divergence free if and only if 

0 = divs(Q;i7) = —d*ai^ = —d*a„ + d*d9 = divs(Q;^) + As6*, (5.1.5) 

where d* is the L^-adjoint of d on S. The above equation is equivalent to: 

As^ = divs(Q;^). (5.1.6) 

This is a Poisson equation. By the divergence theorem, J divs(Q;,y) dAj^ = 0. There¬ 
fore Equation (5.1.6) has smooth solutions. Moreover, any two solutions differ by an 
element in the kernel of As, which consists of constant functions on S. □ 

A minimizer of E can be viewed as a straight out direction of the surface S, as it 
tries to level the surface as much as possible. This motivates the following dehnition: 

Definition 5.1.1 (Uniformly Area Expanding Straight Out Flow). Given a surface 
in a time-oriented spacetime {N^,g). Let T be a global timelike vector field. A 
smooth straight out flow of in N is a smooth normal variation E : S x [0, T] —> N 
such that 

0 C 

YE{x,s) = flufix), divj:fia,,J = 0, (f,Us)dAs = 0, (5.1.7) 

where (3 is chosen such that -^{dAg) = dAg; and s e [0,T], := F(S,s), and Ug is 

a outward spacelike unit normal vector field along with zero divergence. 
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On each surface of in Definition 5.1.1, the outward-spacelike vector Ug with zero 
divergence is only unique up to a hyperbolic rotation with angle 6 (see Proposition 

5.3). The additional condition that (T, VaydAg = 0 defines each Vg uniquely. This 

Js, 

condition is due to J. Jauregui. 

E. Malec, M. Mars and W. Simon have studied such flows in [30], and they have 
shown that given a smooth solution {St} to the straight out flow equation (5.1.7), 

> 0 , 

i.e. that Hawking mass is also monotonically non-decreasing (also see [10]). 

The goal in this chapter is to construct spacetimes that admit smooth solutions 
to straight out flow equation. The idea is similar to the construction in Chapter 4. 
Given U := x M, where Bi is the unit ball in we seek a spacetime metric 

g on U that admit a coordinate chart {t,r,6,(j)}, such that g satishes the following 
four conditions: 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 


,_d_ ^ 

'dr’ 39 


'L A\ 

'dr’ d(j) 


0 ; 

0 ; 


Area form of coordinate sphere St^r satisfies dAt^r = sin OdOdcj)] 


— is straight out, i.e. diV(,„Q;e = 0; 
udr 


(5.1.8) 

(5.1.9) 

(5.1.10) 

(5.1.11) 


where cte,, is the connection one-form associated with Cr- 
The following dehnition is analogous to Dehnition 4.1.2: 

Definition 5.1.2 (Straight Out Flow Coordinate Chart). If a spacetime {N^,g) 
admits a coordinate chart {t, r, d,0} such that the four conditions (5.1.8), (5.1.9), 
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(5.1.10) and (5.1.11) are satisfied, then {t,r,9,(l)} is called a straight out coordinate 
chart, and N is called a spacetime that admits a straight out coordinate chart. 

5.2 Construction of Spacetimes That Admit Uniformly Area Expand¬ 
ing Straight Out Flow Coordinate Charts (work in progress) 


Let U = X M. Given a coordinate chart {t, 

r, 9, 0}, dehne a spacetime metric 

g to be 

t 

r 

9 

<t> 


t 

( 2 

d 

e 



r 

9 ■ = 

d 

u‘^ 

0 

0 

(5.2.1) 

9 

e 

0 

a 

c 



\ f 

0 

c 

bj 



for smooth functions a, b, c, d, e, f, u, v on U. g satisfies conditions (5.1.8) and (5.1.9). 
Choosing two of the three variables a, b and c such that ab — = r^ sin^ 9 satishes 

condition (5.1.10). We show that the fourth condition (5.1.11) is a second order 
elliptic PDE in the metric component d. 

The normal bundle of the coordinate sphere St^r has an orthonormal frame 
{er,en}, where and e„ = jj^- The normal vector n is the same as in 

Equation (4.2.7). We compute the condition on g such that Cr is straight out on 
each S^j.. 

Notation: For the rest of this chapter, the induced metric on St,r is denoted as 
gs; V is the Levi-Civita connection with respect to g] \g\ and are determinants 
of g and gs respectively; and the subscript is omitted from ae,. whenever there is 
no confusion. 
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5.2.1 Computation of The Connection One Form and Its Divergence 


In terms of the frame {d9, df)} for T*St^r, write a = aedO + a^drj), where 

(By dehnition 5.1.1) 


Therefore 



(vi( 

^ d9 


Va ( -y- 
\u or 


/ / \ ^ S n 

= \ —I—V -j—r 

or u or n 


I /=! ^ \ I 

Mllnll \ ^ dr’ Mllnl 


n 


n 


= d 

VA^,n 

ee 


(n is normal) 


I P* 

yT,n 


M||n|| \dt 




|9| 


Mllnii u^\gs\ 


^ \g\ 

-L * 


U 


..If tzMl 
"A tel I 




n^\gs\ 

7.I \ 1/2 


(n is perpendicular to |;, ^, ^) 


(By (4.2.8)) 
(By (4.2.9)) 
(l^l = since \g\ < 0) 


TTl \ V2 


I is) 


(5.2.2) 


By a similar computation, we have: 


I A] - 


Therefore: 


a = — 


u 


d(j) ^ 


77l \ 1/2 




zM 

l^sl 


7.1 \ 1/2 


(5.2.3) 


If (^] 


u 


2 ipr 


7.1 \ 1/2 


• de - 


\gs\ J 


d(j). 


(5.2.4) 
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The divergence of a is defined to be the divergence of its dual vector held, which 
is= where 

= gl^ae + = gf ae + gt^ a ^ (5.2.5) 

Recall that 

Definition 5.2.1. For a vector field X = in local coordinates, its diver¬ 

gence with respect to a metric g is given by: 

Therefore the divergence of a along St^r is 



that is: 

divg^o; = /3 0 + cot 9. (5.2.6) 

5.2.2 Straight Out Flow Coordinate Chart: Big Picture 

We expand out Equation (5.2.6), separating higher order derivatives from lower order 
derivatives. 
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div^gO = + (3^ cot 6 

= (fi'f + 5 '?‘^«<a),0 + (gfc^e + 9s^(^<l>),<P + (fi'f «0 + 5'?«0) cot 9 

= {gfae^g + + gfag^^ + g'g^a^^^) (2nd derivative in g) 

+ (gffif^d + g%a^ + gf^^ae + gt^^a^) + {g^^ag + gl^a cot 9 

(1st derivative in g) 

=: (B) + (il) (5.2.7) 


From (5.2.2) and (5.2.3), ag and are first order in the derivative of the space- 
time metric g. Therefore the divergence equation (5.2.6) is second order in the 
derivative of g. In (5.2.7) above, (B) denotes the second derivative terms, and (/I) 
denotes the hrst derivative terms. 

To compute (B), notice that: 


06 

QS ^9,9 



(plug in (5.2.2) for ag) 




1/2 





7Tl \ V2 


zM\ 

\gs\) 


g^s ' ^"erd + lower order derivatives in g. 




Similarly for the other three terms in (B), the highest order derivative terms are 
the terms containing derivatives of the Christoffel symbols. Therefore 


(B) = (C) + (12), 


(5.2.8) 
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where 


- --0 


1 /- 


l77l\ 1/2 



V I//5I J 

1 

f-isiy- 


11^5| J 

1 ( 

-\a\)''^ 


IssP/" 

1 

(-|9l)‘'" 






■hT^t 


chi 


9^t 




6r.,4i ^ 


«" |gsP/2 


( bT'gj. 0 cPgj, 


P) 


(5.2.9) 


and 






■ (Sl'r', + gl*Tl) 


+ 





(as ^sr + s?r,,) 


(5.2.10) 


Notice that all the second derivatives of g in the divergence equation (5.2.7) he 
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in (D). Note that: 

2 i.99t,r 9rt,9 99r,t) 9 i.96r,r 9rr,9 99r,r) 9 i.999,r 9r9,9 99r,9') 

+ lf'^{994>,r + 9r4>,9 ~ 5'6»r,0)^ 



{9*^\9et,r + 9rt,9) + 9^'^9rr,9 + 9^^9e9,r + 9^'^99cj,,r^ 



c^)(e^r + <^,6») — d{ab 


(?)2uufi + u^{cf — be)a^r + ^^(ce 



^(pr 2 ^9(pt^r 9rt^(l) 9(pT,t) ~t“ ^9(pr,r 9rr,<l) 9(pr,r) 9 ^9(j)9^r 9r9,(p 9cpr,9^ 


9 KQcjxp^r 9r<j)^(j) 9<pr,(l) 

1 


= 2 {9^^i94>t,r + 9rt,cp) + 9^'^9rr,^ + 9^^9<p9,r + 9^^94>^,r) 

= (u'^{ab — c^)(/,r + <^,</.) — d{ab — c^)2uu^^ + u^{cf — be)c^r + u'^{ce — af)b^r'] 

2\9\ ^ ' 


We compute (D) as follows: 

(D) = -;^\(f^{99t,r + 9rt,9) + 9^''9rr,9 + 9^^999,r + 9^^^994,,r\ 

Q 

- 2^t\9^t,r + 9rt,4>) + f'grr,^ + f^9^9,r + f‘^9^p,,r],9 

- 2[t\99t,r + 9rt,9) + ^""9^,9 + t^999,r + t'^99<P,r]^^ 

+ 2^9^^^d4,t,r + 9rt,4>) + 9^^9rr,<j> + 9^^9<j,9,r + 9^‘^ 9 4>(l),r] ,<t> 

= 29^^{^99t,r9 + ^9rt,99 ~ (^9ci,t,r9 ~ — cgg^,.^ — cg^^ g^ + ag^^,.^ + 

2 ^ ^9rr,99 ^9rr,cj>9 ^9rr,9(j> ^9rr,(f>tp) ^9 ^999,r9 ^9(l>9,r9 ^999,r<j> 

+ (^94>9,rcj>) + ~j9^'^{dgg^,^g — cg^^.^g — cgg^.^^ + ag^^^j.^) + {E) 

=: (/) + (//) + {III) + {IV) + {E) (5.2.11) 
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where (E) consists of lower (first) order derivatives in g: 

{E) = ■ {ggt,r + 9rt,e) + 9^,9 ■ 9rr,e + 9^,9'999,r + ’ 994,,r\ 

^ 2 + 9rt,4>) + 9^fi ■ ^rr,0 + 9*fi ' 9^9,r + 9^,0 ' 9<j,cf>,r] 

~ 2^9^,I ’ i99t,r + 9rt,e) + 5'*0' 9rr,e + 5'*0 • 999,r + 9^’^ • 99(l>,r] 

+ 2 ^^0 ' ^ 00 , r ] 


= 2 ^^, 9 ib^,r + hdfi - cf^r - cd^^) + -g^*^{-ce^r - cdfi + af^r + ad^^) 

+ ^9^l{2buufi - 2cuu^^) + + 2auu^^) + “ cc,^) 


+ ^9ll{-ca,r + aCr) + ^g^p{bCr - cb^r) + ^9l^{-CC,r + O&r) 

(see (5.2.11)) 

Therefore the fourth condition (5.1.11) becomes: 


0 = divg^a" = /3j + (3‘^^ + fl^ cot e={B) + {11) 

(see (5.2.7)) 

= {C) + {12) + {11) 

(see (5.2.8)) 

■W + («) + (n) 

(see (5.2.9)) 

= ■{{!) +an + {III)+ {iv) + (E)) 

+ {12) + {11) 


(5.2.12) 


The highest derivatives of 9 (second derivatives) lie in (/) + (//) + {III) + {IV). 
5.2.3 Straight Out Flow Coordinate Chart: Complete Form 
In this subsection, (5.2.12) is computed explicitly from left to right. 
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(!)■■ 


Vj^ {(xh — 

(7^) I_I ^^^96 ^f,r9 C6^j-0 cdj Qfp ~\~ (lf^r(l} ^^,00) 


2|9l 

1 ^ I 

I_I ip^,00 2c(i^^?0 “1“ Ctd^ipcj} “1“ bc^fO cf^^O CC^rcf) “t” ^f,rcl)')' 


(5.2.13) 


Therefore 

1 (-|g|)‘^^ fn_-i-hMl! 

ISsP'" ^ ^ ISsP'" 2 151 

{hd_^Q9 ‘2,cdj Q(p ~\~ Cld^fpfp “h hc_^^Q ^f^rO ^f,r(l)) 



(^(Kee - 2crf,e0 + + {he^re - cf^re 


^^,r<p T ^f^r(p)^ * 


(5.2.14) 


(II)-. 


(II) 


—d[ab — c^) 

m 


{b{u^)fi6 - 2c{u^)^e<i> + «(“^),9i0) 


^|4(6ef=),« - 2cC«=).„* + a(«=),«). 


(5.2.15) 


Therefore 


1 (- 1 ^ 1 )^/^ 


{II) = - 


1 (- 1 ^ 1 )^^^ ,2^ o„9„,2n 


\gs?/‘^ 2 | 5 f| 


zij —(&(m ),ee - 2c{u^)^0^ + a{u^ 


1 / -1 d 


2 Vl^sll^ly 


— {b{u )^ee - 2c{u )fi^ + a{u )^^^). (5.2.16) 


{III)-. 


u^{cf — be) 

{III) ^I_I {bO/ r9 ^^,r(h T CZC ^^). 

2 |^| 


(5.2.17) 
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Therefore 


1 (-1^1)^^^ 1 /{cf-he) 




ISsP'^ 


(III) - 


/- lasP'^ 

cf — be 1 


m 


(bci^r6 — CC^r8 — COj^rtp T Cl'C^rip) 


2 

1 / —1 \ cf — be 


— CC^^rO — C-C^^rcf) ClC^r 


Vl£/s||^| 


ab — (? 


(bct'i'g — CCrpQ — CQjj-(p + OiCy 


(5.2.18) 


{ivy. 


{IV) = 


u^(ce — af) 

m 


{bc j.0 cb ri.Q cc^j.fp “ 1 “ Q/b y.fp). 


(5.2.19) 


Therefore 


1 (-IsD'^^ ,,,n 1 (-|9|)‘'“ u^(ce-af) 


V? IssP'^^ 


(IV)- 


|9s|V^ 

ce — af 1 


2151 


(bc^rO ^^,r9 


2 l9sl=/M-|5l)‘/" 


{bc^rO ^b .j.Q cC j.fp ab j.fp) 


1 / —1 \ ce — af 


2 Vis's![s'!y a6-c2 


{bc^j.Q cb^r^Q cc^7.0 T (xS^y.^). (5.2.20) 


We now compute {E). First we compute the derivative of the inverse metric: 

( 1 ) 


a'l = 


“^ISslV {2uu^e\gs\ + M^|fi'sL 0 )|fi'l - M^lfi'sllfi'Le 


l77|2 


Isl /,e 

2 uu^ 0 \gs\ +2u'^\gs\cot9 u‘^\gs\\g\,e 


|7tI2 


(By (A.2.17)) 


u^\9s\ 

Isl 



+ 2 cot 6 — 
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Similarly 


u^\gs\\ { 2 uu^^\gs\+u^\§s\^ )|^| - 

1^1 \W 


2uu,^\gs\ 15-511 S'1,0 

1^1 \W 


(Ifi'sl = 3111^6') 


m2 

9s\ 

1 1 

( 2 ^ - i 

dl 

’A\ 

1^1 ' 

1 « 

\9 

\) 


Tjtr _ 

y,e — 


-d\gs\ 


{-d,9\gs\ - d\gs\,e)\g\ 

\W 


- (-c^l^sDI^L 


-d, 9 \gs\ - 2d\gs\cot 9 d\gs\\g \,9 

m \9? 


(By (A.2.17)) 


d\9s\ 

(+ 2 cot 9 — 

\9\,e 

\9\ 

\ d 

\9\ 


Similarly 


-d\gs\ 

m 


i-d,^\gs\ - 


\9? 


- i-d\9s\)\9\,<p 


-dj9s\ d\gs\\9\,^ 
m \9? 

-d\gs\ (^0 1 ^ 1 , 0 \ 
1^1 \d \g\) 


/ u^{cf - be)'\ (2mm,^( c/ - be) + u^{cf - be )^9 

V m ^ 

2 mm , 0 ( 0 / - be) + u^{cf - be) ^9 u^icf - be)\g\^e 

m W 


u^{cf - be) (c/ - be )^9 _ 

\g\ \ u cf-be \g\ 


u^{cf - be)\g \^9 
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Similarly 


u^{cf - be)\ _ {2uu^^{cf - be) + u^[c] - be)^^)\g\ - u^{cf - he)\g\^^ 


V 1^1 /,</- 

2 uu^^{cf - be) + u^[cj - be)^ 4 , u^[cf - be)\g\^^ 

m w 

u^jcf - be) / ^ {cf - be)^ 4 , _ 

1^1 \ u cf -be \g\ ) ' 


(4) 


fu^ice - af)\ _ {2uu^0{ce - af) + u‘^{ce - af)^0)\g\ - u‘^{ce - af)\g\^0 

V m ),0 ~ W 

2uu^0{ce — af) + u^{ce — af)^0 vf{ce — af)\g\fi 

i " w " 

u^{ce - af) / ^ (ce - a/),g _ 

1^1 \ u ce-af \g\ ) ' 


Similarly 

-t<i> fu‘^{ce-af)\ {2uu^^{ce-af)+u^{ce-af)^^)\g\-u^{ce-af)\g\^^ 

|9| |9|^ 

2uu^^{ce - af) + u^{ce - af)^^ u^{ce - af)\g\^^ 

m \w 

^ u^jce - af) / ^ (ce - a/),^ _ \g\A 
\g\ \ u ce-af \g\ ) ' 
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Therefore 


{E) = ( 2 — + 2cot6' - ) {be^r + hd^ - cf^r - cd^^) 

2 1^1 \ u \g\ 


+ 2 - 


^u'^\9s\ \g\^^ 


u 


2 —- ^ (-ce^r - cdfi + af^r + ad^^) 


+ 


-du\gs\ f dfi 


1^1 


^ + 2 cot 6 * 
d 


If/I ,6 

1^1 


{hufi - cu^^) + 


-du\gs\ f d^ _ 1^ 
d \g\ 


-cu^e 


+ ClU^fj)) + -- 


lu^{cf-be) f u^e , {cf-be)^e \g\fi 


^ M cf — be 


{ba^r — cc^r) 


^2 


1 u^{cf - be) f {cf - fee ),0 \g\,^ 


1^1 


u ^ cf -be \g\ 


-) { — Ca^r + aCr) 


+ 2- 


lu'^{ce-af) f Ufi {ce-af)fi \g\ 


1^1 


^ M ^ ce-af \g\ 


{bc^r — Cb^r) 


+ 2 - 


lu'^{ce-af) f (ce-a/),<^ \g\,^ 


m 


u ^ ce-af \g\ 


) { — CC^r + ab^r)- 


Therefore 

1 (- 1 ^ 1 )^/^ 

\gs\^/^ 


iE) = - 


1 f -1 \ 


1/2 


vif/.ii^i; 


2 — + 2 cot 6 / 
u 


m 


{be^r + bd^e — cf^r — cd^^) 


+ ( 2 ^ - ) {-ce^r - cdfi + af^r + ad^^) 


/ -1 Y^d 
Vl^5||^|y u 


^ + 2 cot 6 * - ] {bu^e - cu^^) + 


d,(j) |f?|,(/i 

d \g\ 


-cufi + au^^) 


+ 2 


I f -I 


Vl^sll^l 


1/2 


cf — be 
ab — 


u cj — be \g' ' 


\g\,^ 


u 


cf — be 


1 ( -1 

■"nissiisi 


1/2 


ce — af 
ab — c^ 


(— ca^r T ®c^f) 

2 ^ ^ {ce-af),e _ ^ 

u ce-af \g\ 


{bc^r — Cb^r) 


+ ( 2^ + ^ ) (-CC, + abr) 

u ce-af \g 
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That is: 



1 ^1,6 

1^1 


^ + 2 cot 6^ — ] {hufi — cu^cj,) 


^ \9\,^ 

u cf — be 


+ 


cf - he 


\ 9 \ 


ah — 

(— ca^r + ac^f) 


2^ (c/ - 


u 


+ 


cf - he 
ce — af 


{ha^r — cc^r) 


ah — 


rv h \ L I L ^f),<l> \9\,4> \ ( \ h \ 

(he r — ch r) -\- I 2 — +- - — — ,_, I ( — cc T' + db 'p) 

U ce-af \g\ ' 


2 ^ ^ (ce - a/),e _ ^ 
u ce-af \g\ 

(5.2.21) 


Next we compute {12). To do that, we hrst have 




= -^u,9 


= -^u^e 


= -^ue 

I'J 



TTl \ V 2 


zm 

\9s\J 


771 \ 1/2 


-I^l\ 


+ I 2 


11 /- 


771 \ 1/2 


\9s\ ),6i 


77l \ -1/2 


-\9\,9\9s\ + 1^11^51,0 


\gs\J u^2\\gs\J \ \gs\‘^ 

_|^|Xl/2 7 /,._,Xl/2 


|i/sl 


^/MV 

2 u 2 


1 /-MV'V2ii,0 


V l^sl J 

MV'' 

V l^sl ) 

1/zMV'' 

V ISsI ) 


\\9s\\9\ 


1/2 


-\9\,9 2|^||^5|cot0\ 

I//5I \9s? J 

(- 1 ^ 1 , 0 + 2 |^| cot0) 


- J ( ^ ) (- 1 ^ 1 , 0 + 2 |^| cot 0) 
u 2 ' ' 


u 


— + cot 6 — 


\,e 


mj' 


(By (A.2.17)) 
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Similarly we have: 





zM 

|t/sl 


T7l \ 1/2 


+ I- 2 


ttI \ 1/2 


\9s\) 


/ l77l\l/2 


= 1 fzM 


I ~\9\,4>\9s\ + \9\\:9s\^ 

u^^\\9s\) \9s\^ 


V? ■* V \3s\ ) \ \g\ ) V Issi 


1/2 


\9\4 


771 \ 1/2 


2 


1/2 


" V \gs\ J 2^2 

^ f - 

V l^sl ) \ u 2 V 1^1 

l77l \ 1/2 




Ut=!') (-IsL*) 


LfzM 

V l£/s| 



2 m ^0 ^ \g\j(p 

u 2\g\ 


Then we need to compute: 

gffl + g’/Tl = ]^(Z - Z) 


1 1 


2 \9s\\9\ 


v?{ab — c^)(e^r + d,d) — d{ab — (Y)2uufi + u^{cf — be)a^r + u^{ce 


c{u‘^{ab — c^)(/,r + d^Y> ~ d{ab — (Y)2uu^^ + u^Yf — be)c^r + u^{ce — af)br)^ 


1 1 


2 IS'sK&e,^ + bdfi - cf^r - cd^Y “ d\gs\{2buu^g - 2cuu^Y 

+ u^{cf — be){ba^r — cCr) + u^{ce — af){bc^r — c6,r)j 

Vj^ 1 / cf — he 

= TrT=7 i^^,r + bdfi - cf^r - cd^Yj - ibu^ - cu^Y) + —-- cc^ 

2 q V u ab — 


+ 


ce — af 


-^{bc^r - cb^y 


ab — c 


af)cr) 
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Similarly: 


gfrl + 


1 1 


■ ^ — c{y?{ab — (?){e^r + d^) — d{ab — c^)2uufi + u^{cf — be)a^r + u^{ce 


2 IgsM 

■ c^r) + a[u^{ab — c^)(/,r + d^ip) — d{ab — (?)2uu^^ + v?{cf — be)c^r + M^(ce — af)b^ 

= ^ 1 \,-, (u‘^\9s\{-ce^r - cdfi + af^r + ad^^) - d\gs\{-2cuu^e + 2auu^^) 

2 \gs\\g\ ^ 

+ V?{cf — be){ — ca^r + CbC,r) + U^{ce — af){ — CC^r + 

1 / 2(i cf — 6c 

= \ ~ ^dfi + af^r + dd^ij)) — —{—cufi + au^fj)) + — -^{—ca^r + 


+ 


ce — af 


-^^{-CCr + ab^r)y 


ab — c 


Putting them all together, we get: 


(^ 2 ) = ^ 


1 f-\g\Vf‘^^.. 


/ V ) \ u 


“+cot« 


mj 2|j|v' 


( {be^r + bdfi - cf^r - cd^^) 


2d cf — be.^ , ce — af.^ , 

— ipUfi - CU^^) + - CCr) + _ ^2 ~ 


"■x 


1 f-\g\y 1^1,X 1 ^ j ^ 

^ vwj W “ 

Of! n f — hp np — n f \ 

--{-cue + au^) + X- oi-car + ac,) + -- U-cCr + a6,)). 

u ab — ab — / 
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Simplifying {12), we get: 


1 / _1 \ 1/2 
(12) = -- 


VI£/5||^I 



2u 


u 


— + cot 6 — 


|gi,ii 

2151 


{he^r + bdfi — cf^r — cd^^) 


— {bu^g - CU^^) + ^ - %{ba^r - CC^r) + ^“ cb r)'] 

u ab — ab — / 


cf - be 


ce — af 


+ 


+ 


\9\,4> 
u 2\g\ 

cf — be 
ab — (? 


(— — cd^Q + Ojf^r + Ojdfjf) — -(— cu^g + ClU^q 


. , ce — Ojf, , 

i—ca r + aCr) - : -^ ( ~CC r + abj 

’ ’ ab — 


(5.2.22) 


Next we compute {11): 

(^1) = i&(^9 + gl‘^ga<^ + gt^^ae + + (gs^^e + cot 6 

= (gsfi + gt^ + e)ae + + gl^ cot e)a 4 , 

^ ~vf (l^) ^ 

^ ~ff 2 ("1^) ((5's!(9 + gf, 4 > + gf cot 6)Vq^ + {g^^Q + gf^^ + g^^ cot j 


where (/la) := {gf^ + gf^^ + gf cot 9)tI^, and (/16) := (c/J^g + g'^g^^ + g^^ cot 9)tI^. 
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Note that 


(/la) = + gf,p + cot ““ c^)(e,r + d^g) - d{ab - c^)2uu^e 

+ v?{cf — be)a^r + u^{ce — af)c^^ 

= -^'^■^{g^sfi + 5's^0 + 9^s cot 6)d^e + {u^\9s\{g^s,e + 9t^,<j, + 9^s cot 0)e^r 

- 2d\gs\uufi{g^lg + gf^^ + g^^ cot 6) + u‘^{cf - be){gf^g + cot 0 )a,^ 

+ u'^ice - af){gl^g + gf,^ + g^^ cot . 

Similarly: 

(^1^) = (glfi + S'sJ, + S'? cot - c^)(/,r + d^^) - d{ab - c^)2uu^^ 

+ u^{cf — be)c^r + u^{ce — af)b^^ 

= y ^(%!e + 9stf> + S^?cot e)d,^ + ^ (u‘^\9s\{9lfi + 9t4> + ^?cot 0)/,, 

- 2 c/|^s|MM,^(^?g + + 4 ”^ cot 9) + u^{cf - be){gfg + + gf cot 0 )c,^ 

+ u^{ce - af){gfg + g^^^ + cot e)b^^. 
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Plug them back into (/I), we get: 



+ ((^ 1 ^ + + 98 cot e)ufi + ( 4^0 + + gl* cot e)u^^ 

H-^ + 9stt> ■*■ ^)c,r j 

H- ^ 2|^| ((fi's^e + 9s^<p + 5'!^ cot 6')c,r + {9s!e + 9stt> 9^s cot 6')6,,, j |. (5.2.23) 

And 

(5.2.23) = - + 9s^,4, + 9^s cot 0)d^g + + gf^^ + g^^^ cot 0)d^^'^ 

- ^ ((£/& + 9s\ + 9^ cot 0 )n ,0 + [g^^g + g^^^ + g^^ cot d)u^^ 

+ ((^S + 9t4> + cot 0 )a,, + (^^(*0 + g^s'^^ + g^^ cot 0 )c,,) 

+ _ ^{ ((S's^e + 9s^^ + 9^ cot 0 ) 0 ,. + {gl'^g + + g^'^ cot 9)b^r 

It can be computed that 

9sfi + 9t^4, + 9s^ cot 6 = ^(^,0 -bcote- c,^); 
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and 


& + 9s'!’^ + 97 cot 9 = + ccot9 + a,^). 


Using the above, we can simplify (/I) as follows: 

1 / —1 \ 

(/I) = - ( 1 —j-piT ) [{bo — bcotO — c ^)d 0 + (— c e + c cot 6 ^ + a 
^ vlfi'siifi'iy ^ ’ 


1 y -1 

+ - 


2 \\9s\\9\ 


1.^0 \( cf-be ce-af 

{bfi — b cot U — c^(f,) I e^r — ^—« H-;-H-^- -c 


u ab — 


/ ^ d, , cf — be ce — af ^ 

+ (~c ^0 + ccot 9 + a^(j)) j f^r ~ ~af) — {7^’^ ~ab — 


ab — 

(5.2.24) 


Ping (5.2.14), (5.2.16), (5.2.18), (5.2.20), (5.2.21), (5.2.22) and (5.2.24) back into 
the divergence free eqnation (5.2.12), we get: 
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2(-tellsl)''" ■ '((/) + (H) + ("I) + (m 


«2 |gs|3/2 


+ (-^)) + (^2) + (/I) 


(By (5.2.12)) 


d 


= {bd^ge - 2cdfi^ + ad^^^) + {be^re - c/,r-e - ce,r<^ + af^r^) - ^{b{u )^gg - 2c( 


u 




cf — he ce — Qif 

+ (2— + 2 cot 9 — ) (be r + bd g — cf r — cd s) 

\ u \q 


+ ( 2 ^ - {-ce^r - cdfi + af^r + ad^^) - ^ 


^ + 2 cot 6^ — ] {bufi — cu^cf) 


, {cf-be),^ ^ 

u cf -be \g\ 


+ 


cf - be 


ah — 

(— ca^^f + cic^T’) 


2 ^ {cf-be),g _ ^ 

M c/ - fee 1^1 

ce — af 


{ba,r — cc,r) 


+ 


ab — 


2 ^ ^ (ce - af),g _ ^ 
M ce - a/ 1^1 


(fee,, - cb,r) + ( 2^ + ^ ) (-CC, + afe,,) 


u ce-af \g\ 


2u 


u 


— + cot 6 — 


\9\fi 

m 



2d 


be,r + bd,g - cf,r - cd,^) - {bu,e - cu,^) 


u 


+ ^ - ^{ba,r - CC,r) + ^ - ^{bc,r - cfe,^)) - ) ■ ( {-Ce,r - cd,g 


2d , , cf — be ce — af 


ab — ’ ’ ab — 

+ af,r + ad,^)-‘^{-cu,g + au,^) + , „ 

u ab — 

+ {b,g — bcot9 — c,^)d,g + {—c,g + ccotfe' + a,^)d,^) 
+ {b,g -bcotO - c,^) ^e,,. - 2^(i + 

+ {-cg + ceote + a^){fr-2^d+%^^''- ' 


/ \ cc Clf, 1 \\ 

i—ca r + aCr) -i - : - t{~CC r + afe r) 

’ ’ ab — ’ ’ / 


C.r + 


ab — 6 


2 i • 


(5.2.25) 







































Notice that the terms underlined, up to a \gs\ factor, equal the Laplacian of d 
along St,r (see equation (A.2.16) in Appendix A.2.3). We can simplify (5.2.25) as: 


(5.2.25) = \gs\Aggd + {be^re - c/,re - ce^r<f> + 


+ 


d cf — bs 

— (b(^U ^^90 2c('U “1“ Ci(^Vj “h ~ CC^j-Q CQj^np “h 

U QjO c 

- cb^re - CCr<j> + db^rcf,) + ^COt 9 - “ <^f,r “ cd^^) 


^|^(-ce,^ - cd^e + af^r + ad^^) 


2d 

u 


^ , fd,4> 2n. 


+ cot 6^ - 7^ {bu^e - cu^^) + - 

d u 2\g\ ) \ d u 


4 > I5'L</> 


+ 


+ 


+ 


+ 


cf — be 


ab — 

(c/-6e),0 1^1,^ 


(c/-fee),0 1^1,0 


cc, 


cf-be 2\g\ 


) { — ca^r + aCr) 


ce — af 
ab — (? 

{ce-af)^^ 1^1,0 


[ce-af)^Q ^,0 u 

---—— - cot ^ {bcr - cb^, 

ce-af 2\g\ J 


— TUZT ) ( —CCr + Cib^r) 


ce-af 2\g\ 

n r n \ f 2 (i c/ — 

+ (o 0 — 0 cot 6* — c a) e r- u 0 H—-— 

’ ’ \ ’ M ’ ab — 


cf — be ce — af 
ab — c^ 


. ^ „ 2d cf-be ce — af^ 

+ (_c,„ + ccotfl + ( U - -u, + 


- ^ ) (-cn,0 + an,^) 


(5.2.26) 


Definition 5.2.2. We set L 2 {g) = L 2 {d,e, f,u,a,b,c) to be the second derivative 
terms in the straight out flow equation, that is: 

L2{g) ■= \gs\Agsd + {be^ro - c/,r0 - + af^rfl) - ^{b{ifl)^ee - 2c(m^),0^ + 

u 


cf — be 
ab — 


^a^x(j) “1* ^ 2 iP^,T'9 ^b^fQ CC^rcf) “t“ 


(5.2.27) 
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With this defintion, we further simplify (5.2.26) as follows: 


(5.2.26) = L 2 {g) + cot + bd^e - - cd^^) - + bd^e - cf^r - cd^^) 


\9\,^^ , , / j j ^ ‘^d\(d,e 2m,0 I^Lm 

- ^|^(-ce,r - cdfi + at,r + ad,^) - — y— - — + cot 9 - ) {bu^e - 

f d (j, 2u (f) /■ , \ 

+ (t - — - 2M j *■“" 


+ 


cf - be 


(of -be), A _^f(cf- be), 151 ., 


ab — c'^ l\ cf — be 2\g\ 


cf-be 2\g\ 


(- CCl^y> + CIC^T*) 

ce-afv f {ce-af)^e \g\,e 


+ 


ab — c‘^l\ ce — af 2\g\ 


— cot 9 ) {bc^r — cb^r) + 


(ce - a /),0 _ \j^ 
ce-af 2\g\ 


( — CC^r + Cib^r) 


, 2d cf — be^ ce — af ^ 2d, „ 

+ b gCr - b gU g -\ - z --0 gQ r -\ - : - Trb gC r — -\ - bu g COt 9 

’ ’ u ’ ’ ab — c^ ’ ’ ab — c^ ’ ’ u _ 

cf — be ce — af 2d cf — be 

^ ^ba^f cot 9 —bc^Y- cot 9 “1” c^fpu^g ^ 

ce — af 2d cf — be ce — af 

- ~ ^ ,6y + c}Veotd 

2d ^ cf — be ^ ce — af ^ „ „ 2d 

- cu^ cot 9 + — -r-cc y cot 9 + — - -cb^r cot 9 + a ^/,y-a ^ 

u ab — c-^ ab — c^ u 


cf — be ce — af 

H- ; -y®,<*C y H- - - -a ^b r- 

ab — c^ ab — c‘^ 


(5.2.28) 


90 






































( 5 . 2 . 28 ) = L2{g) + cot 9{bd^0 — cd^^) 



{he^r + bdfi - cf^r - cd^^) 



cdfi + af^r + 


2d 

u 



2^ 

u 


+ cot 6 



{bufi - cu^^) + 



^ _ MA. 

u 2\g\) 


{-cufi + au^^) 


cf -be V ( (c/ - be)fi 
ab — I \ cf ~ be 


IsU 

2|9l 



{ba^r 


( {cf - be) 

V cf- 



(— CCL^r “I” CLC^r') 


ce — af 
ab — 


/ (ce - a/),g 
\ ce — af 


laU 

m 



{bc^r 


cb,f) + 


( (ce - af)^^ 
\ ce — af 



{ — CC^r + Cib^r) 


2d 

“t” {b^oc ^p c^fpe^p c^of^p a^fpf^p^ {b^oU Q c^pxi^Q c oCfp “1” a^fpv,^(jfj 

u 

cf — be CC — Qjf 

”t” i w{b^oa^p a^pC^fp e^pC^Q "f* a^pc^p'^ “1” ^ ^{b^oc^p e^pC^fp b^pC^Q 4* a^pb^p^ 


ab — 


2d cf — be ^ 

H-cot (Jibu 0 — cu p) - - - - cot tliba p — cCp 

u ’ ab — 


A cote {be p 
ab-c^ ^ ’ 


cb^p) 

( 5 . 2 . 29 ) 
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Finally the above can be simplified to: 


L 2 {g) + cot 9{bdfi - cd^^) - + bd^ - cf^r 


cd^^) - -^(-ce^r - cd^e + aj\r 


+ CLd^fj)) 


2 d 

u 


d,e ‘^u^e /, N , [d,^ 2m,^ , 

-d- — -m) ‘ ^ It - — - 2^ I 


+ 


cf — he 
ab — 


{cf-be)fi \g\,e „ ^n\(u \g\,^ 

2 cot 9 I [ba^r — CC^r) + 


cf-be 2\g\ 


cf-be 2\g\ 


(—Cflj. + OjCj-^ 


+ 


ce — af 
ab — 


_ |k _ 2 ^te) (tc. - c6.) + ( - M; 

ce-af 2\g\ J \ ce - af 2\g 


-CC^r + Cib^r) 


(b^O^^r r) {b,6^,6 C^QCfj) 

u 


cf — be 


^ ce af 

“t ^,0^,r) —{bgC^^ 


^,r^,0 b^rpC^Q a^fpby^f 

(5.2.30) 


= : \gs\^gsd + F{d,d'). 

Therefore we have the following characterization of the fourth condition (5.1.11): 

Proposition 5.4. The fourth condition (5.1.11) is a second order elliptic PDE in 
d: 

Ag^d + G{d,d') = 0, (5.2.31) 


where G : = 


— = ^.. 
|gs| r* sin^ 9 ' 


If Equation (5.2.31) is solvable, then many spacetimes that admit straight out flow 
coordinate chart exist, even beyond spherically symmetric examples. A necessary 
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condition for (5.2.31) to be solvable if 



G{d, d!) dAt^r 


0 . 


(5.2.32) 


We conjecture that this is always the case. The verihcation of (5.2.32) is still work 
in progress. 
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6 


Conclusions and Open Problems 


We have constructed many examples of non-spherically symmetric, non-static space- 
times that admit smooth global solutions to inverse mean curvature vector flow. 
Prior to our work, such solutions were only known in spherically symmetric and 
static spacetimes. Our work seems to suggest that spacetimes that admit inverse 
mean curvature vector flow solutions might exist generically. However, this more 
general problem is still open: 

Problem 6.1. Given an arbitrary spacetime. Can we always find a “right” initial 
surface such that inverse mean curvature vector flow starting with this surface exists 
for all time? 

Going to the big picture of relating local and global notions of mass, it is still 
unknown that: 

Problem 6.2. Given a spacetime that is sufficiently asymptotically fiat (e.g. 

Schwarzschild outside a compact set). Does the Hawking of inverse mean curva¬ 
ture vector flow surfaces approach the total mass of the spacetime? 

A natural next step is to consider the following problem: 
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Problem 6.3. Given a spherically symmetric spacetime that admits an inverse mean 
curvature vector flow coordinate chart. Consider a perturbation of the spacetime 
metric. Does the perturbed metric admit an inverse mean curvature vector flow 
coordinate chart as well? 

We conjecture that this is always the case for some Minkowski spacetimes: 

Conjecture 6.1. Given Minkowski space with inverse mean curvature vector flow 
coordinate chart that can be smoothly extended to the boundary, consider a pertur¬ 
bation of the spacetime metric. The resulting spacetime still admits inverse mean 
curvature vector flow solutions (in a single spacelike hypersurface) that exist for all 
time. 
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Appendix A 


Geometric Calculations 


A.l Ricci, Scalar and Einstein Curvature of Spherically Symmetric 
Spacetime 


In this subsection, we are going to compute the Ricci, Scalar and Einstein curvature 
of the spherically symmetric space time (iV^,^), with coordinates (t,r,6^,0), such 
that g has the local coordinate representation as in (3.4.1): 


9 = 


t r 9 

—n^(t,r) 0 0 

0 u‘^it,r) 0 


0 
0 

0 0 0 


0 


0 


(A.1.1) 


0 sin^ 6 j 


Note that {dt, C, A, d^] form a local frame of the tangent bundle. We assume that 
the connection on {N^,g) is the Levi-Civita connection. The Einstein summation 
convention will be used, and colons will always denote the coordinate chart derivative 
whereas semicolons will denote covariant derivatives.In the following, we will use 
Latin letters and so on to be indices taking values in {t,r,9,(j)}. Recall the 
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(Christoffel symbols) 


following formulas from Riemannian geometry: 
1 

— ( 

2' 


+ 9ii,j 


9ij,l) 


Ric., - Ric(3., d,) = r‘, - + rLr’” - (Ricd curvature) 

R = tr^Ric (scalar curvature) 

We will write the Christoffel symbols in four matrices r*,r'’,r^, and as in Alan 
Parry’s survey paper [[35]]. Now we compute them in sequence, using the fact that 
the metric g is diagonal in our coordinates. 

1. The computation of Ph 

• r*, = y^\gu,t) = B2vv, = 

• = Kt = l9^\9rt,t + 9tt,r “ 9tr,t) = |^(2nn,^) = 

• ^\e = = l9*\9et,t + 9tt,e - 9te,t) = 0. 

• = P^^ = -^g {g^t,t + 9tt,(j> ~ 9t(j>,t) — 0. 

. P^, = \g^\2grt,r - grr,t) = 1^{2UU,) = 

• Fj.^ = Fg,, = ^9 {.9dt,r + 9rt,6 9r0,t) ~ 0. 

• ^r<p ~ ^(pr ~ 29 i94it,r + 9rt,4> ~ 9r(f),t) = 0. 

• ^69 = \ 9 ^^{ 2 ggtfi — gee,t) = 0. 

• ^00 = F ^0 = \g^\gp,t,e + get^p - 9e<p,t) = 0 . 

• = \9^\‘^94>t,p - gH,t) = 0. 


Therefore we have 


F* = 


t r 9 (j) 

^ VjI Vr 0 0^ 


yu ^ 0 0 

V 


0 0 0 0 
0 0 0 oy 


(A.1.2) 
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2. The computation of T’’: 


• r[t = - gtt,r) = \^{2vv^r) = 

• ^tr = Kt = lf''{9rr,t + 9tr,r “ 9tr,r) = 

• = ^et = \9'’'’{.9er,t + 9tr,e - 9te,r) = 0. 

• = r^t = + 9tr,<P - 9t4>,r) = 0. 

• r^r = l9''''9rr,r = l^{2uU,r) = ^• 

• Ke = ^9r = \9''''{9er,r + 9rTfi “ 9re,r) = 0. 

• = \9'''’{.g4>r,r + grr,^ - 9r4>,r) = 0. 

• ^le = \{‘^9er,e - gee,r) = |^(-2r) = -^• 

• ^64> = = \g''’'{.g4>r,e + ger,^ - ge<t>,r) = o. 

• - 9<P<P,r) = l^i-2rsm^e) = 

Therefore we have 


r sin^ 9 


r" = 


/ 


9 

0 


VV^r 

u 

Q 

U U 

0 0 -T 

0 0 


0 - 


0 

0 

0 

r sin^ 6 


3. The computation of T®: 

. Tf, = y^%2gte,t-gu,o)=0. 

• rt = r^t = \9^%9te,r + gre,t - 9ttfi) = 0. 

• = \g^^{gee,t + gte,e - 9te,e) = 0. 

• = r^t = lg^%g<t>e,t + gte,<p - 9^,9) = 0. 


(A.1.3) 
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• Kr = \g^\‘^gre,r “ = 0. 

• = \g^^{gee,r + grOfi — gre,e) = 

• = lg^^{g^e,r + gre,<t> — gr4,,e) = 0. 

_ p6 _ 1 „d6 _ ri 

• ^ ee - 2 g gss,e - U. 

• ^6<j> = = lg^^ig<t>e,e + gee,!!, — ge<f>,e) = 0. 

• = \g^%2g^e^^ - g<p<p,e) = |^(-2r2 sin6'cos^) = -sinO cos 6. 

Therefore we have 


/ 


= 


9 


t r 9 
0 0 0 
0 0 i 


0 


0 


0 

0 

0 


0 y 0 0 0 — sin 0 cos 9 j 


4. The computation of T'^: 

• rft = lg^^{2gt^^t - gtt,4>) = 0- 

• rt = + gt<i>,T - gtr, 4 ,) = o. 

• ^te = ^et = \g'^'^i.ge<p,t + gt<t>,e — gte,<i,) = 0 . 

• = lg^^igH,t + gt^,<f> - gt^,<f>) = o. 

• rt = \g‘^‘^{ 2 gr 4 >^r - grr,<p) = 0. 

• rt = rt = \g‘^‘^{gd<t>,r + gr<p,6 - gTe,<t>) = o. 

• rt = rt = W^ign.r + grc^,^ - gr4>,<p) = sin^ 9 = ^. 

• rt = |fi''^‘^(2fi'6»0,6» — ge9,4>) = 0. 

• rt = rt = \g^^{gHfi + ge^,^ - ge<^,4>) = l^:2^r^‘2sm9cos9 


(A.1.4) 


cos 6 
sin 0 * 
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• = l9'^'^94>4>,4> = 0 - 


Therefore we have 


Y'f> = 


t r 
0 0 
0 0 
0 0 


9 

0 

0 

0 


COS 6 

sin 6 


(j)\0 i ^ 0 / 

^ \ r sin 0 / 


(A.1.5) 


Now we compute each of the components of the Ricci curvature with respect to 
the basis {5*, 5^, d^]. 


1 . 


Ricti = Ric(5t, dt) 


pr _ pr , pr pt , pr pr- _ pr- pt _ pr pr , p6l pr , p0 pr 
tt,r ■*' tr,t ' rt-*- tt ' ^ rr^ tt tr ^ tr^ tr ' ^ 9r^ tt ' ^ (j>r^ tt 




UtVt 
+ + 
U V 


U^r VV^r 
U 


VV^r V^r 
U? V 


U^t U,t 2 VV^r 

u u r 


1 





2 

H—uUi 
r 


Ur 1 

-+ - 

u 


U,tV,t 

V 



2 . 


RiCir 


Ric(5t, dr) = 



, pt pr _ pt pr , pe pr , y<p pr 

' ^ tr^ tr ^ rr^ tt ' ^ 0r^ tr ' ^ <f>r^ tr 



V J ,r V U 


UU^t VV,r 
^2 


r u r u 


2^ 
r u 


3 . 


RiCie = Ricet = 


pt _ pr _ 

tt,0 ^ tr,0 ~ 


0 
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4. 

Rici,^ = Ric<^t = = o 


5. 


RiCr-r 


r* — r* + r* r* + r* r’’ — r* r* _ r* r'’ _ r® 4 - r® r’’ — r® 

■'■ rr,t ■*■ rt,r ~ ■'■ tt^ ~ -l ^y. -l ± ± r0,r ~ ■*■ rr ■*■ rO^ rO 


■p(/> . p</i pr 

^ r4>^r ' ^ (pr^ rr 


y4> 

^ rcp^ r(p 



^ V,tUU^t 
\ V J ,r V v"^ 


V,r U,r 
V U 


V,r V,r 
V V 


UU^t U,t 2 U^r 

U r U 


V rr 2Ur Vt 1 {UU tt V rU r 

= -^-1-^-^-( -^-1-^^ 

V r u V \ V u 


6 . 


Ric^e = Ric^r = 0 


7. 


Ric^,^ = Ric^r = 0 


8 . 

Ricee 


r* 4- r’’ 4- _ r’’ — r'^ 4 - r’’ — r'^ 

-'- ir-*- 06 ^ ^ 60,r ^ ^ rr^ 00 ^ 06^ 6r ^ 6<p,6 ^ ^ (f)r^ 06 ^ 6(1)^ 0<j) 


r—r (~^\ 'U,r—r —r 1 

) ,r u u"^ r 



1 —r 
r 



2 


1 \ Ur V r 1 
— + r— - r— - 

V 


9. 


Ricg^ = Ric^e = 0 


10 . 


Ric00 


pi pr I pr I pr pr _ pr p</> , p6» i F® _ F^ F*?^ 

^ tr^ <P4> ' ^ (p(p-,'f ' ^ T’r-^ (p(p ^ (fxf)^ (pr ' ^ <P<Pi9 ' ^ (p(p ^ (p(p^ (pO 


-r sin^ 9 


+ 




-r sin^ 9 




u^r-rsm^9 . . „ „ „cos6' 

H- 5 -1- (— sin 9 cos 9) g + sin 9 cos 9 -- 

u ’ sin 9 



sin^ 9 + r sin^ 9 



r sin^ 9 


v,r 1 

V 
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Now we can compute the scalar curvature R of this metric as follows: 

R = tr^Ric = g^^Rictt + g^^RiCrr + g^^Ricoe + gf'^'^Ric^^ {g is diagonal) 


-11/ 2 \ U^r V^tU^t U^tt 

VV^rr + -VV^r - nVV^r H- 

\ r J V u u 



V rr 2Ur V t 1 

-^-1-^-"I- 

V r u V'^ V 


uu^tt 

V 


+ 




V,r 1 

r* ^ 

V 


1 V^rr 2 1 V^r U,r V,r U^t V,t U,tt 1 1 V,rr 2 U^r U,r V,r 

V? V r V? V V u u v"^ v? v r v 


Jl + ^ (1 _ Jl ) + ^ 

u u \ ^2 / ^ ^3 r V 


2 V^rr 
V 


~r Z 

U-^ V 


u 


Utt ^ ^Ur 

2^-^-^ 

I " rj I O 


u 


r w 


A 1 v^r 2 

r V? V 



1 




That is 


R 


2 V rr UrV r UiV t U u ^ r 

-^ + 2 —--2 — —+ 2 —-^- - 

V V u u r 


A 1 v^r 2 
r V? V r2 



. (A.1.6) 


The Einstein curvature tensor G is given by 

G = Ric — ^Rg- (A.1.7) 

In our case of the spacetime metric (A.1.1), we compute the components of G 
with respect to the basis {dt, dr, de, d^]\ 


1 . 


1 1 / 2 

Gtt = G{dt, dt) = Ricti - -Rgtt = ^ vv^rr + -vv^r 




^VV^r 


H— (- U tt j + —R 

u \ V / 2 

2Ur 0 In/, 1 
= 1 -^ 


r w 




(A.1.8) 
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2 . 


f^ic 


-'tr 


1 2ut 

-Rgtr = RlCtr =-^ 

2 r u 


(A.1.9) 


3. 


4. 


Gte — Get — Ricte — ~^Rgte — Ricte — 0 


Gt(j) G (pt RRt(/) 2RR£(/) 0 


5. 


Grr = Ric rr - ^Rg rr 


V^rr 2 U^r V^t 

V r U 


1 fUUtt VrUr\ U 

UU^t + - —^-- - 

V \ V u / 1 


2 


2 v^r vr ^ 1 

rj^ /jj iy>‘ 2 i i^‘ 2 i 


(A.1.10) 


6 . 


G ro G Qr Ric^’^? RiCj’^ 0 


7. 


G y'fp G RiCj-(^ ^RQy-fp Ric^ 0 


Gee = RiC (?0 - -Rgee =1-^ - r — — - —R 




u . 


v^r 1 


U-^ 


V 


2 

r V. 




= ——— r^-A— + G — -R — G-^— — + r\— (A.1.11) 


V 


V 


u 


u 


U'^ 


V 


9 . 


G^^ = Ric,^,;!, - -Rg4>4> = sin^ 9 ■ Ric^e 


-Rr^ sin 9 = sin^ 9 ■ Gee 


= sin^ 9 


/ V rr ^UrV r t V t oU tt ^ Ur l'yr\ 

I-^- r ——^—h r ——^ r — - r ——I- r-^ 

\u'^ V V u u v"^ V J 

(A.1.12) 
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A.2 Calculations in Inverse Mean Curvature Vector Flow Coordinates 

A.2.1 Determinant of the Spaeetime Metric and its Inverse in Inverse Mean Cur¬ 
vature Vector Flow Coordinates 

Given a matrix d, its {i,j)th cofactor is Cij := where Mij is determinant 

of the matrix obtained by deleting the ith row and the jth column of A. The adjoint 
matrix adj(d) is dehned as adj(d)jj := Cji = If A is invertible, then 

= ddf^adj(d). 

Let A be the matrix representation of the spacetime metric g as in (4.2.1). 


adj(d)ii = Cii = Mil = u^{ab — C) (A.2.1) 

adj(d)i 2 = C 21 = -M 21 = -d{ab - C) (A.2.2) 

adj(d)i 3 = C 31 = M 31 = u^{cf - be) (A. 2 .3) 

adj(d)i 4 = C 41 = —M 41 = u^{ce — af) (A. 2 .4) 

adj(d )22 = C 22 = M 22 = -v‘^{ab - C) + /(ce - af) + e(c/ - be) (A.2.5) 

adj(d )23 = C 32 = -M 32 = -d{cf - be) (A. 2 . 6 ) 

adj(d )24 = C 42 = M 42 = —d{ce — af) (A. 2 .7) 

adj(d )33 = C 33 = M 33 = -u^v% - - bd"^ (A. 2 . 8 ) 

adj ( 21)34 = C 43 = —M 43 = + u'^ef + (A. 2 .9) 

adj ( 21)44 = C 44 = M 44 = —u^v'^a — — ad^ (A. 2 . 10 ) 


det(5') = d2lC2l + d22C22 = —dM2i + M22 

= —df{ab — C) + — v‘^{ab — C) + /(ce — af) + e(c/ — 6e)] 

= (—— d^){ab — c^) + eu^{cf — be) + fu‘^{ce — af) (A.2.11) 

These prove Lemma 4.5. 
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A.2.2 Computation of (n, ri) 


Lemma A.l. 


<n, n> = 


det{g) 


det{g) 


\ 9 \ 


u^{ab — C‘) u^det{gs) ' u^\gs\'‘ 


(A.2.12) 


where we set |5f5| := det{gs). 


Proof. We simply compute that 

, . / d —d d cf — be d ce — af d d —d d cf — be d 


+ 


dt dr ab — C dO ab — C dcf dt dr ab — C d9 

ce — af d 
ab — d(j) 


. d d . —d. d d . cf — be. d d . ce — af , d d . 

dP dt dP dr a6 — dP dO ab — c^ dt ’ dp 

-d d d d'^ d d 

dr' dt dr' dr 


+ 


+ 


cf — be d d 


ab — c^ do' dt 
ce — af . d d 


( aTi’ al) + 0 + 


cf — be^^^d ^ [cf — be){ce — af) ^ d 

{ab - c2)2 ^do'dp^ 


ab-c‘^J ^dO'dO^'^ 


ab — c^ dp' dt 


C) + 0 + 


{ab — c^)"^ ^dp'dO^ \ab — c^ 


= -v^ + 


-d^ 


+ 


ab — c‘ 


■{2{cf - be)e + 2{ce - af)f) + 


-d^ d^ 


<- -> 

^dp' dp^ 


^ ^ (ce - af){cf - be)c ^ / cf - be V ^ ^ / ce - q/ V ^ 


{ab — c^y 


ab — (? 


ab — c^ 


= : -v^ + 


-d^ 


+ 


ab — c^ 


{2cef — 2be^ + 2cef — 2ap) + (A) 
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where (A) is 


^ _ ^ (ce - af){cf - be)c ^ / cf - be V ^ ^ ce-a/ y ^ 


{ab — c^Y \ab — J \ab — c^ 


2bc^e‘^ — 2ac^ p + 2abcef + ac^f^ + ab^e^ — 2abcef + bYe^ 


{ab — c^Y 
+ a^bf^ — 2abcef^ 

= — - —(2(Yef — 2abcef — bYe^ — aY p + aYe^ + Ybp\ 

[ab — c"*)"' V / 

= ^ ^ 2)2 (2ce/(c^ - a6) - Y{af + 6e^) + ab[be^ + a/^)^ 

^ — 2cef[ab — Y) + [ap + 6e^)(a6 — 


(a6 — 0^)2 
1 


ab — Y 


[af^ + 6e^ — 2cef). 


Plug (A) back into the above, we get: 

—(Y 1 

<n, n> = H- ^ + — - ^(2ce/ - 26e^ + 2cef - 2a f) 

ab — Y 


+ —- -[af +be -2cef) 

ab — Y 


—r/2 1 

= —H- — H— - :^[‘2cef — 26e^ + 2cef — 2ap + ap + be^ — 2ce/) 


m2 ab — c 


,2,,2 


- Y 1 /O X L 2 

+ ^-^(2ce/ - - af) 




ab — c2 


Y[ab — Y) 


-^ — Y){ab — Y) + Y[2cef — be^ — af)'^ (A.2.13) 


Y\9sy 


(A.2.14) 


by Equation (4.2.2). 


□ 
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A.2.3 Laplacian along St^r 

Note that the laplacian of d on the snrface {St^r,gs) is: 

A,,d = div,,(V,,<;) - ^ Yj L, ^ Y id ((9s'*.)V ISSI 








sin 6 \ 36 
1 


fi -r gs 


sin 0 ) + ^ 
d(j) 


,9 + 


sin 6 


= ^2 0 {i9sp(^,e + gl^d^ee + glpd,<t> + gl^d,^e)r^ sin0 + {g^^d^e + gl^d^^)r‘^ cos 6 

+ igt4>d,9 + gfd^e4> + ds^^d^^ + gs^d^^^)r^ sin 0 + 0 ^ 

= {g^gd^ee + 2g^g d^^e + gt^d^^^) (second order derivative in d) 

+ {g^sfidfi + g^s,ed,<j} + gt^^^dfi + gg'^^d^ff,) + {g^^d^ + ( 75 '^ d^^) cot 9 (A.2.15) 

That is: 

^gsd = {g^s d,ed + "^g^s d,(t>e + gt^d^^^) + {g^g^ + gg^^ + g^g cot 9)d^o 
+ + ^ 5 j. + £/?cot 0 )d^ 

= {{hd^ee - ‘2cd^^g + ad^^^) + {h^ - 26 cot 6 -c^^ + 6 cot e)dfi 
+ (~c 6 » + 2 c cot 6 + — c 

\ee - ‘icd^^e + ad^^^) + {h^ - 6 cot 6 * - c^^)dfi + (-c 0 + c cot 6 * + a. 


|£/ 5 | 


(A.2.16) 


For the second to the last eqnality above, we have used the following computa¬ 


tions: 


\gs\,e = (’"^sim 6')_0 = r ^2 sin 6 ^ cos 6 ^ = 2 r^ sim 0 cot 0 = 2\gg\ cot 0. (A.2.17) 
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and 



(^ ^ 

, ^,e\9s\-h\9s\ 


e \9s? 



26 cot 0); 


Similar for the other derivatives of the inverse of gs- 


±_\ 

\9s\J^ Ifi'sl 

(A.2.18) 


A.3 First Variation of Area 


Let be an embedded closed (compact without boundary) hypersurface in a 

Riemannian manifold {M^, g,^/). Endow S with the induced metric. We consider a 
variation of S as follows: 


F : S X (- 5 ,5) —. M, 5 > 0 , (A.3.1) 

such that for all a; e S* := F(S, t), and t e (—<5, 5), 

^F{x, t) = g{x, t)iy{x, t), (A.3.2) 

where ?7 is a smooth function g e C°°(S x (—6,5)), and z/(a;,t) is the unit outward 
normal vector to S* at {x, t). Therefore the variational vector helds along each surface 
Tjt is -^ = gt ■ Let gt be the induced metric on S*, and let V* be the associated 
Levi-Civita connection. Let dat be the corresponding (n — l)-volume form on 
and At the (n — l)-volume. Let Vt be the n-volume enclosed by S*. We shall refer to 
At as the area of S^, and V the volume, in analogy to the case where are surfaces 
in a 3-dimensional manifold. Let Ib and Ht := tr^JI be the second fundamental 
form and the mean curvature of S with respect to u{x,t) respectively. We first 
compute the variation of dat. Let {U;xi,X 2 , ■ ■ ■ , Xn-i] be a local coordinate chart of 
S, then Tit can be locally parametrized as {xi,X 2 , • • • , Xn-i,t] with each hxed t. Let 
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9t = {gt)ijdx'^dx^ be the local representation of the metric on i,j = 1, 2, • • • , n— 1. 

-zirdat = det(oAdx^ a dx"^ • • • a dx"'~^ 
dt dt^ ’ 

=det( 5 ft) ■ trace ^ a ■ ■ ■ a dx'^~^ (A.3.3) 


2 ^Jdei{gt 

= -Y^ciet(^trace ^ a • • • a dx^~^ 


where equation (A.3.3) follows from the identity 


—det(A) = det(A)trace ^A ) , 


(AAA) 


for any square matrix A with entries functions of t. Now 

d d d d _ d d d _ d 

dt dt ^ dx^ ’ dxt ^ ^ dx^ ’ dxt ^ ^ dx^ ’ ^ ^ 


_ d d d _ d 

5^' df dx^ ^dx^ dxj dt 

/= d , . d — , 

= <V^r/z/,—> + <—, V^r/z/> 

dx^ dx^ dx3 


(V is torsion free) 


= 2T]llt ( ^,^ ) . 
dx^ dx^ 


(A.3.5) 


where the last identity follows from the fact that z/ is a normal to the surfaces. 
Notice that Equation (A.3.5) implies that the first derivative of the metric (along 
the variational vector helds) is given by the second fundamental form. Now plug 
(A.3.5) into (A.3.3): 

= -^/det{gt)tmce f • 2g ■ Ib(^, ^) j dx^ a dx‘^ a ■ ■ ■ a dx"'~^ = Htgdat 

(A.3.6) 

Therefore 

3 c ^ r* 

—At = J ~^dat = j Ht{x)g{x,t) dat{x). (A.3.7) 
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A.4 Second Variation of Area 


Now we compute the first variation of mean curvature Hf, which gives rise to the 
second derivative of area. Recall that Ht = gl\llt)ij in local coordinates, i,j = 
1, 2, • • • , n — 1. Thus 

= + (A.4.1) 

Since 0 = §i{9t9T^) = + 9t{§-t97^), we have §-^gi^ = -gT^{§igt)97^■ Thus 

the first term in the above becomes 


(A.4.2) 

We now compute the derivative of the second fundamental form, 
d d _ d 

/= = d , = d , 

= {V ^ V z/, ——r) + {V ^ z/, V ^ —r) 

3* Bxi dx^ dxi Bt 

3 3 3 

= {V_^ V_^Z/, r) + {(V_5 V_^ — V_5_V_a — Vr_5_ + {V Z/, V _5 r) 

dx'^ 3X^ dx'^ df ’dt-^ 30C^ daz^ dt 

- A"' ij> + <^4’ 4"' 4 ^ < V"- vl> 

- (-Ve..,) . T> + T),, T_> + 


where we have used two lemmas, which will be proved below: 


Lemma A.2. V_^z/ is tangential, i.e., (V_^z/, z/) = 0. 

Bx'^ Bx'<- 


Lemma A.3. V^z/ = — where Vst is the surfaee gradient on Sj 


dt 
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Therefore 


9t 


-AEth - V9t(R{t^, A), /y, + hlllltlP 

ex'- ox^ 

d d 

—Ast?7 — rjg^\R(u, ^:^), —r, z/) + ? 7 ||IIt|P (ambient metric g trace) 
dx^ dx^ 


= -Aj:,g - gRiCg{u, u) + g\\llt\ 


where we used 

Lemma A.4. gl\V^ (Vs,h), = ^^tV- 

dx'^ 

and 


Lemma A.5. gl\V= ||//t|p. 

5cr* dxJ 

Combining above, we get 

- g {RiCg{u,u) + \ \llt\\^) =: Ls,, (A.4.3) 

where Ly,^ is called the stability operator of Sj. The second variation of area is then 
given by: 


= J g{x,t){LYtg){x,t)dat{x) +Ht{x) dat{x) + H^g{x,t)‘^ dat{x). 

{AAA) 

Now we verity the above lemmas. 

Proof of Lemma A.2. Since u is the unit outward normal vector held, we have 

0 = z^> = 2<V^ u, z/>. (A.4.5) 

Thus V^z/ is tangential. Similarly, V^z/ is also tangential. □ 

dx'^ 
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Proof of Lemma A.3. Recall that the gradient of a smooth function rj along the sur¬ 
face Si is dehned as 

Vj:,r] := Vf] - (Vr],u)u, (A.4.6) 

that is, the tangential component of the gradient with respect to the ambient metric. 
For any point p e Si, choose geodesic normal coordinates {f/; ci, 62 , • • • , e^} around p 
such that Cl, 62 , • • • , e„_i span TpSi, and = z/. Since is tangential, it suffices 
to show that (V^z/,ei)(p) = (—VstP, ei)(p), for i = 1,2, • • • ,n — 1. Indeed: 

dt 

(V^z/, CiXp) = -<Z/, V^CiXp) = Vei4 ;Xp) (V is torsion free) 

at at (Jt 

= jy){p) - Ve^Z/Xp) 

= -ei(p)(p) (Vg^z/ is tangential) 

n—1 

= <-Vs^p,eiXp). 


Since p is arbitrary, z/ = — VstP, as desired. 

dt 

Proof of Lemma A.4- First note that gl\V_^CVY:tV): 

axt 

But 


□ 

sFVS (VE.rri, 5 fj>. 


sFvy (Ve,.,), X> - i’',.(V(VE,.,)) = div,,(VE.i() - A,:,,,, (A,4,7) 

where V*(VstP) is the covariant derivative of the vector held Vs^p, hence is a (1,1)- 
tensor held. □ 

Proof of Lemma A.5. Dehne vector helds X{i) := V_5_z/ and Y{j) := V^z/. Using 


dx'^ 


dxJ 


n—1 Q n—1 Q 

local coordinates, we can also write X{i) = 

fc=l Z=1 


dx’-' 
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Then 


I III,I r - (yu), 

- 9?(ft)„«A(i)"y(i)'’ 

-gf(X(i),y(])) 

dx3 

as desired. □ 


A. 5 Transformation Formulae of Ricci and Scalar Metric under Con¬ 
formal Change of Metrics 

Given Riemannian manifold {M^,g), recall that the Riemann curvature operator R 
acts on vector helds X,Y,Z e r(TM) as follows: 


R(X, Y)Z = DxDyZ - DyDxZ - D[x,y]Z. 


(A.5.1) 


We can thus define a (1, 3) tensor field R = Rfijdx^ (x) dx^ (x) dx^ (x) such that 




' dx^ ’ dx^ dx^ dx^ 


One can verify that 


RC- = 




It I -pK pm p/c pm 

.• J- A yy. O 7 -1 m -*- ■? 7 • 




(A.5.2) 


Using the metric g, we can define a new (0,4) tensor field Rm := Rkiijdx^®dxh 
dx* (X) dx^ with components Rkuj := gkmRuj- One verihes that 


Rkiij = (R{di,dj)di,dk}. 
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Rm is called the Riemannian curvature tensor with respect to the metric g. The 
Ricci curvature in the direction X e TpM is given by 


Ric(VX):=Afi(V,T)^.X>. 


(A.5.3) 


Therefore the Ricci curvature tensor is given by the (2,4)-contraction of the Riemann 
curvature tensor, i.e., 

Ric., - R.c(T. . g>^RjUk = /'R,,,, - -4, 


' da;* ’ dx’^ dx^ ’ dx^' 


dT^. 

^ ik , -p/c p/ _ pfc p/ 

ki ij ji ik- 


(A.5.4) 


The scalar curvature is 

R := ^f^-’RiCij. (A.5.5) 

Proposition A.6. Let (M**, g) he a Riemannian manifold of dimension n, and p > 0 
is a smooth function on M. Consider the new metric g := pg, i.e., g is a conformal 
change of g. Then the corresponding Ricci curvatures changes in the following way: 


Rij ^ij ' 


n 


(log p) (log p) (log p) J - ]^gij 


Tl — 2 

^^iVlogpl^ 


(A.5.6) 


Proof. We will verity equation (A.5.6) at the center p e M of a geodesic normal 
coordinate neighborhood. Then at p, we have gij = Sij and = 0, and therefore 

gij,k = dk(di, dj) = <r^dm, dj) + (di, Cf^dn) = 0 at p. (A.5.7) 

By (A.5.4), the Ricci curvature of g is given by: 




_ pfc 

^ ik,j 


I -ps pfc _ pfc pfc 
^ ij^ sk ^ ik^ sj' 


(A.5.8) 
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Thus at p, we have: 


Rij = + ^Sikilogp)jk + ^^jfc(logp),ifc - ^gij^kg^\\ogp)^i - 


-gijg^\\ogp)^ik 


^ 1 1 1 
r^fcj + 2^ik{^^Sp),kj + ^hkihgp),ij - -gik,jg'‘\\ogp),i 


+ (^r|i + ^<^i 41 ogp),j + ^(5j,(logp),i - ^gijg^\logp),^ 
^sk + ^^sfc(logp),fc + ^4fc(logp),s - ^PsfcP'''(logp),Z 

r|fc + ^^«(logp),fc + ^Sks(}ogp),i - ^Pifc/'(logp),z j X 


+ ^^«fc(logp),i + ^<5jfc(logp),, - ^gsjg^^(logp)^i 


— Rij + 


2 — n 


(logp)^ij - ^gij(logp)^kk + (product terms), 


(A.5.9) 
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where the red terms cancel, green terms cancel and the bine terms vanish at p. 


(prodnct terms) = -6jJkk(}ogp)^i(}ogp)^s - -4s<5jfc(logp),i(logp),^ 

1 Tt 

+ PjfcP^^(logp),z(logp),i + ^Pijp"'(logp),,(logp),/ 

Tl — 1 1 — Ti 1 

= ^^(logp).i(logp)j + ^^Pp(logp),/(logp),« + ^P*i<^/fc^ifc(logp),z(logp),z 

- |(5iz(5j7(logp),z(logp),z 

Tl — 1 Tl — 1 1 

= ^^(logp).i(logp)j-^Pp(logp),/(logp),« + ^P*i(logp),fc(logp),fc 

- |(logp),i(logp)j 

Tl — 2 2 — Tl 

= ^^(logp).i(logp)j + ^^Pp(logp),/(logp),« (A.5.10) 

Where the red terms cancel. Therefore 

~ 2 — Tl 1 Tl — 2 

Rij = Rij + {log p)^ij - -pp(logp),fcfc + ^^(logp),i(logp)j 

2 — Tl 

+ -^9ij(}ogp)^i(logp)^i 

= Rij - ^^^(logP),ii + ^^^(logP),*(logP)j - \ (^^ff((logp)) + ^^^|V(logp)|J^ 

(A.5.11) 

□ 


Corollary A. 7. Continue from Proposition (A.6), then the scalar curvature ofg is 
given by: 


R 



n — 1 

p2 


\P 


{n 


1)(^ 

4p3 



(A.5.12) 
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Proof. We again compute at p, which is the center of a geodesic normal coordinate 
neighborhood. The result then follows from a simple calculation: 


n — 2 


n — 2 


R := g^^Rij = ( Rij -^(logp),ij + —^(logp),i(logp), 


^(Ag(logp) + ^^^|V(logp)|2) 


1 ^ n — 2 ^ , n — 2 •., ,, , n 

= + ^^^^(logp),i(logp)j - — 


n — 2. 


Ag(logp) + iVlogpl 


= -R- 
P 


n — 2 n 


n ~ 2, 


2p 2p ' ^ 

1 2(n-l) 


+ 7^ ^—iVlogpP - 


= -R- 
P 


2p 


-Ag(logp) + 


4p 

-rR + 3n — 2 
4p 


n{n — 2) 
4p 


Vlogpl 


Vlogpl 


= -R- 
P 


n — 1 f 1 


. I + 
P \P J 


—rr + 2)11 — 2 ( 1 


4p 


P^ 


)P,iP, 


1 n — 1 
= -R - 


P 


= -R 


P P 
1 


P \ P^ 

AgP + 


1 

P 


-rR + 3n — 2, 


+ +P,ii ) + ^ \^p\g 

1„ n —1. 4n — 4 — + 3n — 2 

■|Vp| 


4p3 




P P 
as desired. 


4p3 


□ 


We want to eliminate the gradient term, and for that purpose we need to distin¬ 
guish the following two cases. 

Corollary A.8. Suppose the manifold is of dimension n = 2. Let p := where 
u > t) is a smooth function on M, then 

R = {R - 2Agu). (A.5.13) 
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Proof. We plug p = into equation (A.5.12), and again, we will be computing at 
p e M in a geodesic normal coordinate neighborhood: 

R = p-^R-p-^AgP + p-^\Vp\^ 

= e-^^R - 

= e-^^R - {e^^2u^i2u,i + e^^2u,ii) + 

= e-2“ {R - 2Agu ), 

where the red terms cancel. □ 

4 

Corollary A.9. Now suppose manifold is of dimension n ^ 3, we set p := , for 

u a positive function on M. Then equation (A.5.12) becomes: 

R = u-^ (^Ru - . (A.5.14) 


Proof. The result follows from plugging p = u "-2 into equation (A.5.12), again we 
will be computing at p: 
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R = u "- 2 /?—(n —1 )m "- 2 Ag(M"- 2 )— 




-u "- 2 |V(m"- 2 )|^ 


^ ^ . 4 6 —n 

= u "- 2 /? — [n — l)u "-2 I -^■U"- 2 m^ 


(n- l)(n-6) / 4 

- U "-2 


n — 2 


2(6-n) 

U "-2 U^iU^i 


= u "-2/?—(n —1 )m "-2 


4 /6 — n 8^ 


n — 2 \n — 2 


U "-2 + M"-2M 


,22 


(n —l)(n —6) / 4 ~i 2 +i 2 - 2 » 


n — 2 


u "-2 M 


,2 ‘-^,2 


_4_ „ , ,4 _I1±2 

= M "-2i<—(n —1)--M "-^AgU 

Th ^ 


|2 ( 4 , (n-l)(n-6)/ 4 


-|V«| (n-1) 


n — 2n — 2 


-u "-2 + 


n — 2 


u "-2 


_^4_ ^ 4(n — 1) n + 2 

= u "- 2 i? . 


(n - 2) 


M "-2 Agti — |Vm|^ ^- 


-4(n — l)(n — 6) ^ (n — l)(n — 6)4^ 


(n- 2)2 4(n- 2)2 


_ 2n 

U "-2 


4 ^ 4(n — 1) _r^ ^ 

= U "- 2 /?-^- —U "-^AgU 


_ n + 2 

= u "-2 Ru — 


(n - 2) 

4(n — 1) 


n — 2 


AgU , 


(A.5.15) 


which is desired. 


□ 
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